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APRESENTAÇÃO 


Os estudos multidisciplinares nas ciências exatas e da terra e 
engenharias são primordiais para o desenvolvimento de soluções inovadoras 
para os desafios do mundo atual. A complexidade desses desafios exige 
que profissionais de diferentes áreas executem atividades em conjunto para 
compartilhar seus conhecimentos e perspectivas. Por exemplo, algumas áreas 
com estudos multidisciplinares: (1) Engenharia de software: conhecimentos de 
matemática, computação e engenharia para desenvolver softwares que atendam 
às necessidades de usuários; (2) Geofísica: estuda os fenômenos sísmicos, 
elétricos, gravitacionais e magnéticos, incidentes sobre o planeta e tem várias 
conexões com Geociências. (3) Biofísica: aplica física às pesquisas biológicas, 
estudando os fenômenos da vida desde a escala molecular até sistemas 
complexos; dentre outras áreas fantásticas. 

Ainda, os estudos multidisciplinares oferecem n benefícios, tais como: 
(i) maximizar a compreensão dos problemas complexos: isto é, ao trabalhar 
em equipe, os profissionais de diferentes áreas podem compartilhar seus 
conhecimentos e perspectivas, o que implica uma compreensão mais detalhada 
dos problemas complexos; (ii) desenvolvimento de soluções inovadoras: a 
combinação de diferentes perspectivas técnico-científicas convergem para 
o desenvolvimento de soluções inovadoras; (iii) otimização da comunicação 
e colaboração: ou seja, trabalhar no ambiente multidisciplinar exige que os 
profissionais desenvolvam habilidades de comunicação e colaboração, eficazes. 

Portanto, para difusão e aplicação dos estudos multidisciplinares, as 
universidades e empresas devem investir nos programas e iniciativas de 
fomento a colaboração entre profissionais de diferentes áreas. É necessário 
que as universidades ofereçam os cursos e programas interdisciplinares e as 
empresas utilizem as equipes multidisciplinares para trabalhar em projetos. 
É fácil notar que os estudos multidisciplinares são tendências crescentes nas 
ciências exatas e da terra e engenharias. De fato, eles são essenciais para o 
desenvolvimento de soluções inovadoras para os desafios nos negócios e na 
vida. Diante disso, o livro apresenta uma fundamentação teórico-prática nos 
resultados obtidos pelos diversos autores e coautores no desenvolvimento de 
cada capítulo com conhecimento técnico-científico, autêntico. Além disso, a 
Atena Editora garante qualidade na divulgação científica e excelência, fatores 
essenciais para conquistar o destaque entre as melhores editoras do Brasil. 


Fabrício Moraes de Almeida 
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ABSTRACT - The objective of the research: 
to identify the significance of epidemiological, 
socioeconomic and coverage indicators 
of health services in Rondônia, applying 
the multivariate statistical technique. The 
method: multivariate statistical technique, 
using Factor Analysis and Principal 
Component Analysis method. A sample 
with 121 variables was used, divided 
into two scenarios: Scenario 1, with 121 
epidemiological, socioeconomic and 
health service coverage variables in the 
State of Rondônia; and Scenario 2, with 
42 epidemiological, socioeconomic and 
health service coverage variables with an 
explanation greater than 95%. Results: The 
statistical inferences show the following 
findings: factor (1) called epidemiological 
indicators, presented 69 variables with 
explanation above 70%, factor (2) called 
socioeconomic indicators highlignted 4 
variables with explanation above 70% and 
factor (3) indicators of coverage of health 
services highlighted only 2 (two) variables 
with explanation above 70%. Scenario 2, 
composed of 42 variables with explanation 
above 95% in the application of Factor 
Analysis and Principal Component Analysis 
obtained 34 variables with explanation 
above 70%. Of these, 25 are of the factor 
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(1) epidemiological indicators and presented explanation above 95%. The variable with the 
greatest explanatory power is the first with 99.66% (referring to the % of children aged 10 to 14 
years, with more than 1 year of school delay). The correlation between factors (1), (2) and (3) 
is explained by 66.53% for epidemiological indicators, 18.72% for socioeconomic indicators 
and 7.53% for service coverage indicators. From this study it is possible to infer the almost 
absolute predominance of epidemiological variables and, with less emphasis, of socioeconomic 
variables, which measure access to elementary education in the State of Rondônia. The 
variables of coverage of health services did not present significance. Conclusions: the study 
showed that of the 39 indicators agreed between the State Department of Health of Rondônia 
and the municipalities, 33 (thirty-three) are included among the indicators with statistical 
significance, which validates the choice of method. 

KEYWORDS: Statistical and Health; Indicators; Multivariate Analysis; Rondônia; Brazil 


11 INTRODUCTION 


In the process of building the Brazilian Unified Health System (SUS), develop 
research projects that have in their core the interest of researching political alternatives 
to foster the allocative policy of financial resources of the SUS, with a view to mitigating 
regional inequalities of access to health services, configures as essential, mainly to raise 
to the national manager the need to adopt financing policies that focus on the health needs 
and living conditions of the population, making prevail the maximum constitutional postulate 
of the Constitution of Brazil [1] “health a right of all and duty of the State”. 

Given the proportions of the SUS system, public health policies have demonstrated 
the capacity for innovation with the implementation of the Basic Care Floor, despite the 
fact that the historical and cultural behavior of resource allocation based on production 
prevails, maintaining a process of distributive inequities. We seek to correct this situation 
through instruments such as the Integrated Agreed Programming (PP), implemented since 
2000, which uses a historical series of production to allocate resources in medium and high 
complexity actions, which comprise the SUS. 

This new model of health care is governed by principles that seek to guarantee 
universal, integral, egalitarian, equitable and free access to health, consisting of a network 
organized in a regionalized and hierarchical way, which advocates a single command at each 
level of government, based on a strategy of administrative and operational decentralization 
of health actions and services and their legal provisions constitute the main instruments of 
regulation of the system for the transfer of resources to the health area between the three 
levels of government [2], [3], [4], [5], [6], [7], [8]. 

Considering the magnitude of the Brazilian State, the allocation of financial resources 
in health without observing regional inequalities, especially in the North and Northeast 
Regions, which are doubly peripheral, significantly accentuates the regional inequalities 
that reside in different levels of income concentration, low demographic density and low 
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living conditions, pointed out through health indicators [9]. Thus, health inequalities are 
expressed in the different possibilities of life expectancy at birth, how to live, the ways of 
getting sick and dying for the same class in contrast to the rates of other regions of the 
country [9]. Itis necessary to consider the regional differences regarding the characteristics 
of the installed service network, access to information and strategic planning instruments 
for the better development of health actions and services, and the investment capacity 
of the different States and regions in technological and managerial innovations, among 
other decisive factors in the performance of the critical areas evaluated in the process of 
reviewing financial ceilings [10]. 

It is evident the need to develop studies that point to public financing strategies 
in health that share with the particularities of the peripheral states of the country, and 
can from then on serve as a basis for (re)discussing the allocation of resources in the 
SUS, respecting the geographical, populational, socioeconomic and cultural aspects of 
each region of the country. The peculiar characteristics of the geographical space and the 
demographic distribution, as well as the inevitable problems of access and locomotion, 
to meet differentiated demands it is necessary to foresee modes of remuneration in line 
with the degrees of difficulties and time required in the execution of the work. A typical 
example to demonstrate the importance of this argument is in the vaccination action, 
essential for all Brazilians: for the sake of justice, it cannot be remunerated in the same 
amount and per capita criteria in the Northern region as in the other regions of the country, 
especially if compared with the South, Southeast and Midwest regions of Brazil. That is, we 
must consider the existence of a permanent challenge for managers at the three levels of 
government, in seeking alternatives to conduct the SUS guidelines, in the perspective of 
reducing regional inequalities, which imply in the distribution/allocation of federal financial 
resources that reach the differences and state and regional specificities in the management 
of the system. These challenges must be overcome, concomitantly with the exercise of the 
capacity to identify and solve the main health problems that affect the local population. And 
that it can modify the perspective of results, observing equal treatment for equal needs, 
equal access and equal health. In this sense, Escoda [11] considers the difficulties of the 
SUS as inherent to a social process and points out its cultural, political and technological 
dimensions. Highlighting that due to its social nature, it is a process still under construction. 
In this process of construction, there are so many laws, norms and normative instructions, 
where controls are still being implemented, and where there is a natural lack of service 
providers in the health area, there is also a lack of more consistent information and reliable 
indicators that adequately reflect reality. The correct situational diagnosis can be extremely 
decisive for the population of the State of Rondônia, and may mean the difference between 
access to treatment (life) and complete abandonment, resulting in deaths. 

What can be deduced is that there is a multiplicity of factors involved: its human 
significance, the volume of resources, its social impact and its uniqueness. The manager 
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of national public health policies has focused on strategic lines of discussion for the 
construction of a management pact, pact for life and pact in defense of the SUS, in this 
context the strongest point discussed within the collegiate instances is the reduction of 
regional inequities. Given the breadth and complexity of the problem, the present research 
was limited in an exploratory study of epidemiological, socioeconomic and coverage 
indicators of health services in the State of Rondônia. 

To this end, we sought to apply the multivariate statistical technique, using Factor 
Analysis and the Principal Component Analysis method, with emphasis on identifying the 
significance of the indicators. According to Mingoti [12] Factor Analysis has as its main 
objective to describe the original variability of the random vector X, in terms of a smaller 
number m of random variables, called common factors “and that are related to the original 
vector X through a linear model”. In this model, part of the variability of Xi attributed to 
common factors, with the remainder of the variability of, Y attributed to variables that were 
not included in the model, that is, to random error. In general, what is expected is that the 
original variables Xj, i=1,2,...,p are grouped into subsets of new variables that are mutually 
uncorrelated, and the factor analysis would aim to find these clustering factors. In cases 
where there is a large number of variables measured and correlated with each other, it 
would be possible, from the factor analysis, to identify a smaller number of new alternative 
variables, uncorrelated and that somehow summarize the main information of the original 
variables. These new alternative variables are called factors or latent variables. On the 
other hand, Principal Component Analysis, from the moment the factors are identified, their 
numerical values, called scores, can be obtained for each sample element. Consequently, 
these scores can be used in other analyses involving other statistical techniques, such as 
regression analysis or analysis of variance [12]. Although Factor Analysis can be applied 
to the original variables contained in the X vector, to facilitate understanding we prefer to 
introduce the main concepts of this technique using the original variables Xi, standardized 
by the respective mean and standard deviation. 

Given this reasonable argumentation, the objective of the research is to identify 
the significance of the epidemiological, socioeconomic and coverage indicators of health 
services of the 52 municipalities of the State of Rondônia, applying Factor Analysis and 
Principal Component Analysis. What is expected with the use of these two statistical tools 
can guide future situational diagnoses in health, in order to effectively guide the allocation 
and use of resources. 


21 MATERIALS AND METHOD 


2.1 Materials 


The object of study in question is the epidemiological, socioeconomic and coverage 
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indicators of health services, with emphasis on a contribution to the planning of actions based 
on situational diagnosis in health. The choice of indicators was based on the notes of Mingoti 
[12] to ensure the quality of the sample data “most multivariate statistical techniques use only 
complete observations, that is, if for a sample element, if the value of any variable has been 
lost, it is eliminated from the analysis process”. The 52 (fifty-two) municipalities of the State 
of Rondônia were considered, described in alphabetical order and acronym: Alta Floresta 
D'Oeste (ALFL); Alto Alegre dos Parecis (ALALG); Alto Paraíso (ALPA); Alvorada D'Oeste 
(ALV); Ariquemes (ARQUEMES); Buritis (BUR); Cabixi (CAB); Cacaulândia (CAUC); Cacoal 
(CACOAL); Campo Novo de Rondônia (CNRO); Candeias do Jamari (CJA); Castanheiras 
(CAST); Cerejeiras (CERJ); Chupinguaia (CHUP); Colorado D'Oeste (COLOR); Corumbiara 
(COR); Costa Marques (CMARQ); Cujubim (CUJB); Espigão D'Oeste (ESPIG); Gov. Jorge 
Teixeira (GOVJTEIX); Guajará-Mirim (GMRIM); Itapuã D'Oeste (ITAPUÃ); Jaru (JARÚ); Ji- 
Paraná (JI-PR); Machadinho D'Oeste (MACH); Ministro Andreazza (MANDREAZ); Mirante 
da Serra (MSERRA); Monte Negro (MNEGRO); Nova Brasilândia (NBRA); Nova Mamoré 
(NMA); Nova União (NU); Novo Horizonte (NHO); Ouro Preto D'Oeste (OPRETO); Parecis 
(PARECIS); Pimenta Bueno (PBUENO); Pimenteiras D'Oeste (PIMEN); Porto Velho (PVH); 
Presidente Médici (PMÉD); Primavera de Rondônia (PRIMARO); Rio Crespo (RCRESPO); 
Rolim de Moura (RLM); Santa Luzia D'Oeste (SLUZIA); São Felipe D'Oeste (SFELIPE); 
São Francisco do Guaporé (SFCO); São Miguel do Guaporé (SMIGUEL); Seringueiras 
(SERING); Teixeiropólis (TEIX); Theobroma (THOB); Urupá (URUPÁ); Vale do Anarí 
(VANARÍ); Vale do Paraíso (VPARAÍSO); Vilhena (VILHENA). 

In order to establish fidelity, the data collection was concentrated in the main 
databases considered as official in the operationalization of the Health System at the national 
level, being: National Registry of Health Establishments (CNES), Department of Informatics 
of SUS (Datasus), Integrated Agreed Programming (PPI), Integrated Health Information 
Network (RIPSA), Public Health Budget Information System (SIOPS), Brazilian Institute 
of Geography and Statistics (IBGE) and United Nations Development Program (UNDP), 
considering the information present between the period of 6 (six) years. The construction 
of the database took place through isolated capture in each information system through 
online access via ADSL with connectivity and home accessibility. Each variable presented 
was extracted from the isolated and main system, and grouped in a table of the Microsoft 
Excel program. 

The data capture period was approximately one (1) year due to the large number 
of variables. Once the process of capturing and exploring the information systems was 
exhausted, only the variable that presented continuous information in the 52 (fifty-two) 
municipalities of Rondônia was attributed as useful. A total of 57 socioeconomic variables, 
54 epidemiological variables and 10 variables of health service coverage were collected from 
the databases of the 52 municipalities of the state of Rondônia. The computer instrument of 
choice for the treatment of data statistically was the software STATISTICA version (10) for 
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presenting analytical capacity of considerable expression. 
2.2 Method 


2.2.1 Description of the search procedures 


The first procedure to arrive at the statistical method was to carefully observe the 
data collected in order to identify distortions in the data collected, thus ensuring continuity of 
the information contained in the variables in the 52 municipalities. The research in question 
is classified into quantitative and qualitative in view of the presence of information in the 
form of data and the representative components of each variable, but it is an exploratory 
research in the sense that the intention is only to explore the statistical amplitude of each 
variable in the scenario State of Rondônia. 

Data were statistically treated using the Factor Analysis technique and the Principal 
Component Analysis method. The procedure at first was to apply the technique and the 
method in the 121 grouped variables and after applying only in the group of variables 
with explanation greater than 95%, divided into two scenarios: Scenario 1: 121 variables; 
epidemiological, socioeconomic and coverage of health services in the State of Rondônia; 
and Scenario 2: 42 variables; epidemiological, socioeconomic and coverage of health 
services with an explanation greater than 95% in the State of Rondônia. 


2.2.2 Factor Analysis 


Mingoti [12] proposes: 


2.2.2.1 Factor Analysis Model via Correlation Matrix 


Be Koxi a random vector with vector of means y, where H=(U,M, «4 P)', covariance 
matrix E, and correlation matrix Eai Are Z, [ (X- 4) / O]the original standardized 
variables, where yu, and o,, represent respectively the mean and standard deviation of 
the variable X, i = 1,2...p. In this case, the matrix se is the matrix of covariances of the 
random vector Z = (2, Z,...Z,)”. 

The factor analysis model constructed from the theoretical correlation matrix Pa it 
is a model that linearly relates the standardized variables and the common factors that, at 
first, are unknown. The equations of the model are given by: 


LM Pt. HF tE, 
Lola Ft Lo Pot... dm FF Es 
5 : (4.1) 


Lo =1 F,+ Io Ft. + Di Fr E, 
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LEAL F,t HF +O 


1=ii 12 2 im m i 


LF 


de Ft Lo Foto ta FA rO, 


In matrix notation, the model (4.1) can be expressed by: 


D(X-u)=LF+ & (4.2) 
where: 
DG - tax, - mall IE Lelleofa 
X — upa = Epa = Fma = Loxm = 
ly lo ce hm 
a ta o los 
Edi » 
lo 0.0 
0 1/0. 0 0 
Dow = ; 
0 o boss tio, 


In this model, based on Mingoti [12], F 


mx1 


is a random vector containing m factors, 
also called latent variables, that describe the elements of the population under study and 
are not observable, 1<m=p, that is, they cannot be measured a priori. Therefore, the factor 
analysis model assumes that the Z, variables are linearly related to new random variables 
F, j=1,2,..., m, that will need to be identified in some way. The vector Epa is a vector 
of random errors and corresponds to measurement errors and Z, variation, which is not 
explained by ordinary factors F, j = 1,2,...m, included in the template. The coefficient of o 
commontly called loading, is the coefficient of the i-th standardized variable Z, in the j-th 
factor Fj and represents the degree of linear relationship between Z, and F, j=1,2,...m. 


2.2.2.2 Orthogonal factor model 


Some assumptions are necessary to operationalize the estimation of the model in 
(4.1). Let's assume that: 


() E[F,,.,] = 0, which implies that E[F] = 0, j = 1,2,..., m, that is, all factors have an 


average equal to zero: 
1 0 «= O 


o 1 : 0 
(ii) Var LP quai = ii = . l - 
o 0 1 
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that is, all F, factors are uncorrelated factors and have variances equal to 1, 


(iii) E [£,,,] = 0, which implies that E [£] = 0, j = 1,2,...,p, that is, all errors have 
averages equal to zero; yo 0 0 
1 
(iv) Var [e ,,] = 'Y pp = , V. o 
0 0 cp, 


that is, Var [E] =y, and Cov (e, E)= 0, “i x j which is to say that the errors are 
uncorrelated with each other and do not necessarily have the same variance; 
Fm)=E(£F)=0 
Assumption (v) implies that the vectors £ and F represent two distinct sources of 


(v) The vectors £,,, and F,,, are independent. So, Cov (£ ,, 
variation, related to the standardized variables Z,, there is no relationship between these 
sources of information. A factorial model with assumptions (i)-(v) It is called orthogonal, 
where orthogonality refers to the fact that the m factors are orthogonal to each other. 

An immediate consequence of assumptions (i)-(v) is related to the structure of the 
theoretical correlation matrix Poa When the orthogonal model is assumed, the matrix as 


can be reparameterized in the form: 
Ea =LL+wW (4.3) 
This comes from the fact that: 


Px = Var(Z) = Var (LF + €) 
Fon Var(LF) + Var(e) = LIL + py = LL +y 

where lis the dimension identity matrix pxp. 

The objective of factor analysis according to Wíchern; Johnson [13] is to find the 
matrices Loem and a that can represent the matrix ia for a given value of m, less than 
the number of original variables p. Unfortunately, there are many correlation matrices Pas 
which cannot be broken down into form LL + W'for a value of m much less than p. 


In (4.4) it is possible to better visualize the format of the matrices involved in the 
decomposition given in (4.3) of the correlation matrix. 
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Poxp = , ; a E BE ga 4 (4.4) 


m m 


Aly L, 
j= 


The implications of decomposition (4.3) are presented below: 
GOvarlZ = [+04 ity= hi +4/,, where hi =1atot..+n,i= 1,2..., Pp. 


which means that the variance of Z,, it is decomposed into two parts. The first, 


2 
Lo; 


ja 


denoted by hã is the variability of Z, explicada pelos m fatores incluídos no modelo fatorial. 
This part of the variability is called “commonality”, a name given in function of the fact 
that the factors F, j = 1, 2, .... m appear in all model equations (4.1) and variables Z, 
have a common source of variation. The second, denoted by Wy, is the variability part of Z, 
associated only with random error e, which is specific to each variable Z,. This part of the 
variability is called “uniqueness” or “specific variance.” As the variables Z, have variances 
equal to 1, it follows that hi+ y = 1. 


GI) col 7) tt o ik= 1,2, pi ék 

(1) Cov(Z, F) = Lom and therefore, Cov (ZF) = Corr (ZF) = | o ear 0 

= 1;2,.;m. 

This comes from the fact that 

Cov (Z, F) = Cov (LF+e, F) = Cov(LF, F)+ Cov (e,F) = Cov(LF, F) =L 

Thus, one can use the L , matrix in the search for understanding and interpretation 
of the factors Foj= 1,2 assi 

To operationalize factor analysis in practice, we must first use mechanisms to 
estimate the value of m. From the estimated value of m we can then estimate the matrices 
Lui and Vox" 

(i,) Regarding the total variance, the proportion explained by the factor F, is given by: 


PVTE = ÉS * (4.5) 
p 
and the most representative factors in the model are those with higher values of 
(4.5). It is common to express the values in (4.5) in percentage. 
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4.2.2.3 Estimation of the Number of Factors 


The first step in conducting factor analysis is to estimate the theoretical correlation 
matrix Em through the sample correlation matrix Roxo as was done in principal component 
analysis. For the estimation of m, itis enough to extract the eigenvalues from the (5,48 matrix 
and sort them in descending order.It is observed, then, which eigenvalues are the most 
important in terms of numerical quantity, using the following criteria: 

Criterion 1: the analysis of the proportion of total variance related to each eigenvalue 
A, given by A/p, i = 1, 2, ..., p. Those eigenvalues that represent higher proportions of the 
total variance remain, and therefore the value of m will be equal to the number of eigenvalues 
retained; Criterion 2: the comparison of the numerical value of À, with the value 1,i = 1,2, 
«= P. The value of m will be equal to the number of eigenvalues À, greater than or equalto 1. 
The basic idea of this criterion is to maintain in the system new dimensions that represent at 
least the variance information of an original variable. This criterion was proposed by Kaiser 
(1958); Criterion 3: observation of Cattell's scree-plot [14], which shows the values of l; 
ordered in descending order. By this criterion, the graph looks for a “jump point”, which will 
be representing a decrease in importance in relation to the total variance. The value of m 
would then be equal to the number of eigenvalues prior to the “jump point”. This criterion is 
equivalent to Criterion 1. 

Suppose, for example, that we had p=6 and the eigenvalues À, =2,24, À,=1,38, A,= 
1,21, À,=0,63, À,=0,41, À,=0,13. 

In this case, by Criterion 2, m would be estimated to be equal to 3. The same 
suggestion would be indicated by the scree-plot shown in Figure 1. 

The criteria described take into account only the numerical magnitude of the 
eigenvalues. An adequate choice of the value of m must, however, take into account the 
interpretability of the factors and the principle of parsimony, i.e. the description of the 
variability structure of the random vector Z with a small number of factors. 

It is important to emphasize that the orthogonal factorial model should only be 
applied in situations in which the original variables are correlated with each other, because, 
otherwise, each factor will be related to only one original variable, making the value of m 


equal to p. 
2 
ER 
a T T "E T T 
1 4 3 a 5 [=] 
ordem do fator 
Figure 1: Scree-plot. 
Source: Mingoti [12]. 
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2.2.2.4 Methods of Estimation of Matrices Lom and Vox 


Once the value of mis chosen, it is possible to estimate the matrices Loem e Wo 
Initially, the principal components method will be introduced, commonly used as an 
exploratory analysis of the data, in terms of the underlying factors, because it does not 
require information or assumptions about the probability distribution of the random vector Z. 
Subsequently, the principal factors method will be presented, which is a refinement of the 
principal components method, and the maximum likelihood method, which is indicated only 
when the random vector Z has multivariate normal distribution. 

2.2.2.4.1 Principal component method for estimating matrices L,,, and 

The principal components method works as follows: for each eigenvalue À,, i = 1,2..., 
mretained in the estimation of the value of m, as discussed in section 3, is the corresponding 
normalized eigenvector e where e= The matrices Land W,, 

shall be estimated respectively by: 

Im!" (4.6) 

Idiag (4.7) 

where diag (.) denotes the diagonal matrix. Thus, the iss matrix has the main 
diagonal equal to the elements of the main diagonal of the matrix “ 

The basic idea of this procedure lies in the application of the spectral decomposition 
theoremtothe R,, matrix. By this theorem, the sample correlation matrix can be decomposed 
as a sum of p matrices, each related to an eigenvalue of the matrix R,,,. For a fixed m value 
one has that: 


Ro, =Thea= Sites Saco (4.8) 
Thus, an approximation to the matrix LL* will be given by: 
To build the matrix !., one can consider using the following matrix: 


x et= R pxp (4.9) 
fome] 


Since the matrix in (4.9) is not diagonal, it cannot be used completely for the 
estimation of Vox However, one can consider its diagonal. Thus, the matrix of specific 
variances is estimated as given in (4.7). 

Considering this form of estimation, the original sample correlation matrix Fa is 
being approximated by: 


R ty (4.10) 
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and the residual matrix from the adjustment of the factorial model will be given by: 

MRES =R-“” (4.11) 

The residual matrix can serve as a criterion for evaluating the quality of fit of the 
factorial model. Ideally, your values should be close to zero. However, this matrix is only null 
when the value of m is equal to p, which in practice is not the desired solution. 


2.2.3 Principal factors method for estimating matrices L 


From Mingoti [12] it is perceived that another method can be used for the estimation 
of the matrices of loadings and specific variances introduced by Thompson (1934) and is 
called the method of principal factors or iterative principal components. For it to be used, it 
is necessary that the value of m has already been estimated by some criterion. The basic 
idea is to proceed to a refinement of the estimates of L,., and Wp,,, generated by the method 
of principal components. Considering the model P = LL'+ Wy, where P is the theoretical 


correlation matrix of the random vector of interest X. Then we have: 


É Pa Pe 
Du hi Dos Pis 


DO, DO,» Ds si h, 


onde h2= 1 -w, i=1,2...., p são as comunalidades. Suponha que se estime a 
matriz LL por R* dada por: 


A Fo No CC Fo 
Fa he -Fa o Pa 

Rpo= | o PR mes (4.13) 
Para fia di) 


p 
where tm ich) are initial estimates of commonalities Wishi- hj). Using the 
principal components method, one has that: 


L = (4.14) 


where are eigenvalues of R*e are the respective normalized einenvectors. From 
the matrix L”, we have, therefore, new estimates of commonalities Wisin), which are 
then placed on the main diagonal of the matrix in (4.13), and the procedure of estimating 
the matrix L*, by the principal components method, is repeated again. The algorithm is 
proceeded until such time as the differences between the commonalities of two successive 
interactions are negligible. Problems may occur during the execution of the algorithm, 
making it difficult for the iterative process to converge. 
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For example, eventually at some stage of the procedure, some eigenvalues of the R* 
matrix may be negative, as these depend on the initial estimates of the commonalities. This 
creates an inconsistency with the fact that one is estimating a defined positive matrix, as is 
the case with the correlation matrix. 

There may, also according to Rencher [15] be the Heywood problem, occasioned by 
the fact that at some stage of the iterative process some estimate “ti can be greater than 
1, which generates a negative estimate of * à , and thus inconsistent with the definition of 


variance. 


2.2.3.1 Principal Component Analysis 


The use of Principal Component Analysis can under the developed conditions be 
understood as a process of analysis of statistical data, since it allows the methodological 
“walk” from beginning to end, that is, it allows the researcher to generate the expectation of 
results. Starting from the need to evidence the technique and method applied to the process 
of analysis of statistical data, and in attention to the objectives of the research, in order 
to leave it as a proposal for a real contribution to the Health System in Rondônia - more 
precisely to the planning area. From the task of explaining the method applied follows the 
demonstration of the statistical mode of how to obtain statistically results in a given sample. 

Mingoti [12] proposes: 


Be X=(X, X, ... Xo)'a random vector with measure vector yu = (4, M,... u,)'and matrix 
of covariances 2 ,,. 

Are A 24,2... A the eigenvalues of the array Exp with their standard eigenvectors 
€, Oy 3 E, that is, the eigenvectors e, satisfy the following conditions: 

() e/e=0forallizj; 

(i) ee =1forali=1,2,...,p; 

(iii) Zwme=A4 epforalli=1,2,...,p 

where the eigenvector e, is denoted by e, = (e,e,,...; €,)”. Consider the random vector 
Y=0"'X, where Da is the orthogonal matrix of dimension pxp, consisting of the normalized 


eigenvectors of the matrix Z,,, that is, 


e, es: RR OR e, 
Ci E» +: e,» 

Omw= [0 0 CCC “|-leere, (3.1) 
e, Op cc Cpm 


The vector Y is composed of p linear combinations of the random variables of the 


vector X, has vector means equal to and O” y and matrix of covariances A... which is a 


pxp? 
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diagonal matrix, whose elements are equal to a ,=A, ,= 1, 2, ..., pthatis. 


Therefore, the random variables that constitute the Y vector are uncorrelated with 
each other. Thus, the idea of using linear combinations in Y arises as an alternative way 
of representing the covariance structure of the vector X, trying to obtain a reduction of the 
space of variables, going from the dimension p, to a dimension K less than p. Therefore, 
instead of using the original random vector in the data analysis, using the K main linear 
combinations. The random vectors, X and Y, have the same total variance and the same 
generalized variance, and the Ysector has the advantage of being composed of uncorrelated 
random variables, thus facilitating the joint interpretation of these. The following are some 
important definitions. 


Definition 1: The j-th principal component of the matrix 2x = 1,2,...,p is defined as: 
Y=e X=e, X,+6,X, +... +6,X, (3.2) 


The hope and variance of the component Y, are, respectively, equal to: 
ElY] = e, = CM, tCpH,+... +e,H, 
Var [Y] =e, “Zpxp e = A, 


being CovlY, Y]=0,jxk. Each eigenvalue Aj represents the variance of a principal 
component Y, Since the eigenvalues are ordered in descending order, the first component 
is the one with the greatest variability and the p-th is the one with the lowest. 

Definition 2: The proportion of the total variance of X that is explained by the j-th 
principal component is defined as: 


var ly | 1 A 


e = 
Variância Total de X Traço (7) px Ai 
ia 


(3.3) 


By the spectral composition theorem, the total and generalized variances of the 
random vector X can be described through the total variance and generalized variance of 
the random vector Y since, 


p p 
traço (Zpxp) = BX-A =>. À, , where, =Var [Xi = 2 ssP 
iã iã 
Thus, in terms of these two global measures of variation, the vectors X and Y are 
equivalent. In general, the ratio in (3.3) is multiplied by 100, indicating the result as a 
percentage. It is evident that the first principal component has the highest proportion of 
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explanation of the total variance of X. 
Definition 3: The proportion of the total variance that is explained by the first k 
compon 


k 
Full dg? 
Variância Total de X “Toi, ) ) À; 


(3.4) 


lf the first k principal components explain a large part of the total variance of the 
vector X, it may be restricted the focus of attention only to the random vector (Y, Y, ... 
Y)”. In this way, a k-dimensional set of random variables can be examined, rather than a 
p-dimensional set, without losing much information about the original structure of variances 
and covariances of the X vector. 

By the spectral decomposition theorem, by restricting the focus of attention only to 
the first k principal components, the matrix of covariances E, pis being approximated by the 
formula in (3.5): 


k 
Das no 2; e;e; (3.5) 


Each portion of the sum in (3.5) involves a matrix of dimension pxp corresponding 
only to the information of the j-th principal component, j = 1,2,...,k Thus, the original 
variability system of the vector X is being approximated by the sum of k matrices, each 
representing the variability system related to a component. When k=p, one has that the 
covariance matrix Ex p IS accurately reproduced by the sum of matrices rationed to the 
principal components, tati is, 


p 
Eai; 
o 


Definition 4: Another way to define the principal components is presented below. 
Consider the following system of linear combinations of X consisting of p equations of the 
type: 

Y=a,X=a,X +ra,X+...+a,X,i= 1,2. 

Thus, one has to: 

Var (Y) = à',Z,xp à; 

Cov(Y, Y)=a, E papi j=1,2..sp 

Suppose one wants to find the values of the coefficients aij such that a”, a,= 1, so that 
the linear combinations Y,, Y,, ..., Va were uncorrelated with each other and had maximum 
variance. 

So it can be demonstrated by Wichern; Johnson [13] that the maximum variance of Y, 
=”, X, under the restriction a”, a, = 1, is equal to À, and is obtained when a, =e,, that is, the 
normalized eigenvector corresponding to À,, being Y, called the first principal component. 
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The maximum variance of Y, = a”, X, under the restrictions a”, a, = 1 and cov(Y, Y,)=0, is 
equal to A, and is obtained when a, = e, the normalized eigenvector corresponding to À, 
being Y, called the second main component. The maximum variance of Y, = a”, X, under 
the restrictions a',a, = 1 and cov(Y, Y,=cov (Y, Y,) = 0, is equal to À, being Y, called 
the third main component. Following the procedure, the maximum variance of Y, = a”, X, 
under the restrictions a”,a,=1 e cov(Y, Y,)=0 j<i, is obtained when a, =, that is, the 
normalized eigenvector corresponding to the eigenvalue A, being Y, called the i-th principal 
component. Thus, the principal components p are constructed, being unique, except for the 
signal exchange of all their coefficients. So, for example, if Yi is a principal component, --Y, 
will also be a principal component of order i. 

Due to the very form of construction, the first principal component is always the most 
representative in terms of total variance and the p-th is always the least representative. 
Figure 2 shows a graphical illustration of the principal components in the case of two 
variables. As can be seen in the said figure, each point in the coordinate system Y, and Y,. 


Y2 


X1 


Figure 2: Geometric interpretation of principal component analysis — p = 2. 


Source: Mingoti [12]. 


2.2.3.2 Principal Component Estimation: Covariance Matrix 


In practice, the matrix E xpiS unknown and needs to be estimated through the sample 
data collected In general, the matrix Empis estimated by the matrix of sample covariances 
pxp" 
Let À, A, «.. A, be the eigenvalues of the Spxp matrix, and let ê,, É, ..., Ê, be the 
corresponding normalized eigenvectors. Then, the estimated j-th principal component is 
defined by: 


V=6/X=8,X,+8,X, +... 48, X 5 j=1,200p (3.6) 
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Some properties of the sample principal components are presented below.. 

Property 1: The estimated variance of Y, is equal to A, sf =d,225D 

Property 2: The covariance between the components Y, and Y, is equal to zero, for 
all k x j, which means that these components are uncorrelated. 


Property 3: The total variance explained by the j-th sample component is given by: 
A; - A; 
rats EA 


ia 


Property 4: The estimated correlation between the j-th main sample component and 
the random variable X, i= 1, 2,...,p is given by: 


= ent; 
v Su 
where S, is the sample variance of the random variable X,. 
Property 5: By the spectral decomposition theorem, the matrix of covariances S,,, 
can be expressed as: 


=V).0.0; 
S5/;6;8; 


or approximated by (3.7) if only the first k sample principal components are used. 


So Die) (3.7) 


In practical terms, to make use of the k main sample components considered most 
relevant in data analysis, it is necessary to calculate their numerical values for each sample 


element, values called component scores. 


31 RESULTS AND DISCUSSION 


3.1 SCENARIO 1 = 121 EPIDEMIOLOGICAL, SOCIOECONOMIC AND HEALTH 
SERVICE COVERAGE VARIABLES IN THE STATE OF RONDONIA 

In the application of the Factor Analysis technique and the Principal Component 
Analysis method, the findings combine with the theoretical foundations and demonstrate 
the importance of developing scientific studies in contribution to the management of the 
health system and in contribution to the decision-making process at the time of distribution 
of technical resources, financial and human. 

In the application of the Factor Analysis technique and the Principal Component 
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Analysis method in the variables, with 121 variables grouped together. The statistical 
inferences show the following findings: factor (1), with an explanation greater than 70%, 
epidemiological indicators presented 46 variables, socioeconomic indicators 17 variables 
and indicators of health service coverage 6 variables. Therefore, with the application of the 
Factor Analysis technique and the Principal Component Analysis method in the variables, 
with 121 grouped variables, 69 variables with an explanation greater than 70% were found. 

Factor (2) socioeconomic indicators highlighted 4 variables extracted from the 
application of Factor Analysis and Principal Component Analysis: Percentage of children 
aged 10 to 14 years with more than 1st year of school delay; Percentage of people in 
households with piped water; Percentage of children aged 10 to 14 years with less than 
4 years of schooling; Percentage of children aged 7 to 14 out of school. The factor (3) 
indicators of coverage of health services highlighted only 2 (two) variables extracted 
from the application of Factor Analysis and Principal Component Analysis, higher than 
70%: Percentage of government income; Percentage of people with more than 50% of 
government income. 

Statistically in the context of factorial exploration, the factor (1) epidemiological 
indicators presented a variance explanation of 69.90%, the factor (2) socioeconomic 
indicators of 16.47% and factor (3) indicators of coverage of health services 6.46%. The 
set of factors accumulated eigenvalues of 92.8% of explanation extracted from the Principal 
Component Analysis. Considering what Mingoti [16] points out “the eigenvalues are ordered 
in descending order, the first component is the one with the greatest variability and the very 
bad one is the one with the lowest”, the first main component is revealed in the study as 
the epidemiological indicators, the second the socioeconomic indicators and the third the 
indicators of coverage of health services. 

Considering also the notes of Mingoti [12] the first main component is placed in the 
multivariate analysis in this study as of greater representativeness, and by inserted context 
of greater importance in view of the object of study. Table 1 shows the statistical findings. 


; % Total umulative umulative 
Values Eigenvalue Varian Einenvalu % 
À á E 69,90479 57,77255 69,90479 57,7725 
ao ds o 16,47150 13,61281 86,37628 71,3853 
FiSpnieam niçes 6,46530 5,34322 92,84158 76,72858 


Coverage Indicators 


Table 1: Demonstration of the eigenvalues in the application of PCA in the 121 variables in the year of 
study 


Source: Prepared from secondary data. 


The eigenvalues of correlation can also be explained through the total variance that 
has the ability to synthesize the global variance of the multivariate distribution, since this is 
the sum of the variances of all the variables involved in the X vector. 
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Another finding that draws attention is the coordination of factors with cases, which 
for this study were analyzed 121 variables and 52 cases, demonstrates that, considering 
the 52 cases, there are in factor (1) 13 cases with negative inferences, in factor (2) 27 cases 
and in factor (3) 24 cases. Itis worth mentioning the health regions hosted by the following 
municipalities: Porto Velho, Ariquemes, Ji-Paraná, Cacoal, Rolim de Moura and Vilhena, 
emphasizing that only Porto Velho (capital) presents itself negatively in a decreasing 
situation in the three factors. 


Caos F (1) Epidemiological F (2) Socioeconomic F (3) Health Services 
Indicators Indicators Coverage Indicators 
5 ARQUEMES -7,8181 4,2101 2,9247178 
9 CACOAL -8,8542 5,4915 2,88012 
24 JI-PR -11,4610 4,6260 2,01097 
37 PVH -47,8678 -11,4336 -2,23120 
41 RLM -4,9203 6,8001 0,65169 
52 VILHENA -9,0195 8,5311 1,33950 


Table 2: Demonstration of negative and positive cases (city of health region) based on correlations, 
extracted from Principal Component Analysis, with 121 variables. 


Source: Prepared from secondary data. 


3.2 SCENARIO 2 = 42 EPIDEMIOLOGICAL, SOCIOECONOMIC AND HEALTH 
SERVICE COVERAGE VARIABLES WITH AN EXPLANATION ABOVE 95% IN 
THE STATE OF RONDÔNIA 

From the exploration of the 121 variables, it was necessary to consider the notes of 
Hair et al [16], who propose, in the use of factor analysis, to adopt criteria of percentage of 
variance in order to obtain practical significance for the determined factors, considering it 
desirable that the level of explanation of variance be 95%. 

According to Hair et al [16], if in the interpretation of factor analysis it is inferred that 
the: 

1) factorial matrix: works as an aid in the process of choosing the number of 
factors, this non-rotated matrix, demonstrates the particular combination of original variables 
explaining more the variation in the data as a whole than any other linear combination of 
variables. Therefore, the first factor can be seen as the best summary of linear relationships 
displayed in the data; 

1) the second factor is orthogonal to the first: (...) thus, the second factor can 
be defined as the linear combination of variables that explains most of the residual variance 
after the effect of the first factor has been removed from the data; 

2) Factor loadings: is the correlation of each variable with each factor. The 
loads indicate the degree of correspondence between the variable and the factor, and larger 
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loads make the variable representative of the factor. 

Scenario 2 represents the findings of the application of Factor Analysis and Principal 
Component Analysis in variables with explanation above 95% in the State of Rondônia. 
Scenario 2 is extracted from Scenario 1, in order to obtain greater clarity of data exploration 
and to envision a greater possibility of contributing to the situational diagnosis in health in 
the planning area. 

The variables with an explanation greater than 95 % are shown in Table 3, as follows: 
They were found 25 Socioeconomic Variables, 8 Epidemiological Variables and 3 Health 
Services Coverage Variables. 


25 Socioeconomic Variables 


% children aged 10 to 14 years older than 1 year | % children aged 10 to 14 Ra with less than 4 
behind in school (CUMATE 5) Ye ars of schooling ( (CQUATES) 

% of children 7-14 years with 1 year of E aged 7 to 14 years in primary school 
delay in school (CG 7-14C1AATRASADO 

% children Ras holds with per capita domá EN child dada 7-14 years out of school (CÇ7- 


E nro of government Re] 

income (PRMGOV) % of government income (RENDAGOV) 

% of peóple in households with electricity and % of E in households with piped water 
RA PDENERG) (GDAGENC) 

en 25 years of age and older omen 15 years of age and older 
[MULAC2SANOS lo to 14 years (POP1O- a 
14ANOS) 9 y Population up to 1 year of age (POP1ANO) 

beplation 25 years of age or older Population 15 years of age and older 
Eanes a 65 Desdrs of age or older a a ; 

(POPAC6SANOS) 5-year-old population (POP<5ANOS) 

Urban population (POPTOTAL) Total population (POPRURAL) 

Probability of survival up to 60 years Probability of survival up to 40 years 

Dna pal Human Devolbpinéit Index Longevity nina Human Development Index (IDHM) 
IDHM-LON Life expectancy at birth (ESPVIDAQNA 

8 Epidemiological Variables 

Hospitalization infectious-parasitic diseases Hospitalization for digestive diseases 
DOENINFEC-PARASIT RGE A 

Hospitalizations for pregnancy complications Injuries from external causes LESCAUEXTER 
COMPGRAVIDEZ . 
General death OBITOGERAL Infant death OB-INF 


Death from External Causes OB-CAUS-EXTER Infant death by residence OB-INF-RES 


3 Health Services Coverage Variables 
Examination of the cervix EXPREVENTIVO Te neRa Ear una anro mn 


NVC/ZCONSPN 
Medical visit number MÉDIAVISDOM 


Table 3: Demonstration of socioeconomic variables with explanation greater than 95% in the application 
of Factor Analysis and Principal Component Analysis for the State of Rondônia. 


Source: Prepared from secondary data. 


E 
promo dida 
cigecono 


Table 4: Demonstration of the concentration of variables by factor and percentage of 
explanation. 


Source: Prepared from secondary data. 
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As for the relationship between the first and the last variable, according to Mingoti 
[12] “variance serves to measure the degree of linear relationship between two variables”. 
Thus, in the application of PCA in the 42 variables above 95%, the variable with the greatest 
explanatory power is the first with 99.66% (% children aged 10 to 14 years with more than 
1 year of school delay; CUMATES), from the second and successively until the 8th, 99.99% 
of explanation is obtained, as shown in Table 5. 

he variability associated with the random error E, which is specific to each variable, 


can also be highlighted. 


Variables Eigenvalue da o no Gu 
1 CUMATES 1,661123E+11 99,66815 1,661123E+11 99,6681 
2 CQUATES 5,354738E +08 0,32129 1,666478E+11 99,9894 
3 CENFUN 1,149565E+07 0,00690 1,666593E+11 99,9963 
4 CDVIQMS 3,275427E+06 0,00197 1,666626E+11 99,9983 
5 RENDAGOV 1,776124E+06 0,00107 1,666643E+11 99,9994 
6 PRMGOV 4,844556E+05 0,00029 1,666648E+11 99,9997 
7 PDAGENC 3,033109E+05 0,00018 1,666651E+11 99,9998 
8 PDENERG 1,187135E+05 0,00007 1,666653E+11 99,9999 


Table 5: Demonstration of covariance eigenvalues of matrices (factors) 


Source: Prepared from secondary data. 


There is a correlation between factors (1), (2) and (3), which explains 66.53% for 
epidemiological indicators, 18.72% for socioeconomic indicators and 7.53% for service 
coverage indicators. The finding corroborates what Hair et al [16] say: “the first factor can 
be seen as the best summary of linear relationships displayed in the data”. 

According to Reis [17] and Mingoti [12], commonality is “the total amount of variance 
that an original variable shares with all the other variables included in the analysis” or “the 
variables have a source of variation in common”. In Scenario (1) the set of 121 variables 
(epidemiological, socioeconomic and coverage of health services) presented 69 variables 
with commonalities above 70%. In Scenario (2), the set of 42 variables (epidemiological, 
socioeconomic and coverage of health services with explanation above 95%) presented 16 
variables with commonalities above 95%. It is noteworthy that the variables of Scenario (2) 
are included in scenario (1). 
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F(1) 


varios | Epidomiotoical GlQlndkatos Servicos Mute, 
NMAC15ANOS 0,983867 0,996511 0,996612 1,000000 
MULAC25ANOS 0,985501 0,996481 0,996542 1,000000 
POP1ANO 0,965979 0,994573 0,995047 0,999977 
POP10-14ANOS 0,978427 0,996382 0,996495 0,999994 
POPAC15ANOS 0,982552 0,995959 0,996183 1,000000 
POPAC25ANOS 0,984372 0,995766 0,995977 1,000000 
POP<5ANOS 0,974340 0,995863 0,996087 0,999988 
POPAC65ANOS 0,979914 0,980371 0,980497 0,999872 


Table 6: Demonstration of commonalities in variables above 95% explanation 


Source: Prepared from secondary data. 


41 CONCLUSIONS 


There is an almost absolute predominance of epidemiological variables in factor (1), 
which represents an explanation above 70%. This fact is external, thus the affinity of the 
object of the work in the face of being represented by the health problems, especially the 
causes of hospitalizations and deaths. 

Regarding factor (2) of socioeconomic variables, there was a predominance of 04 
variables that tend to express the quality of the education system of Rondônia in view of the 
presence of indicators that measure access to elementary education. On the other hand, 
factor (3) - of the variables of coverage of health services - did not present significance for 
the study. 

For the 52 municipalities in question, the situation of the Municipality of Porto 
Velho (Capital of the State of Rondônia) stands out, which presents a negative correlation 
of (-47.86%), which can be justified by the presence of diseases common to the other 
municipalities, for example, dengue, tuberculosis, leprosy, injuries due to external causes, 
among others, here observing only numerical values. Only the municipality of Porto Velho 
(Capital) contributes with 72.84% in Factor (1) and 17.64% in Factor (2). This fact can be 
explained by the characteristic of the city in concentrating the services of medium and 
high complexity in health and being references for all the municipalities of the State in the 
following services: urgency and emergency in adults (Hospital Estadual Pronto Socorro 
João Paulo — Il), major surgical treatment (Hospital de Base Dr. Ari Pinheiro), pediatric 
urgency and emergency (Cosme e Damião Children's Hospital) and in infectious-contagious 
diseases (Hospital Centro de Medicina Tropical de Rondônia). In this case, some variables 
such as hospitalization and death can directly influence the contribution of the municipality. 

Scenario (2) aimed to stratify the result of Scenario (1) in order to ensure better 
clarity to the study. In the application of factor analysis and principal component analysis 
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in 42 variables with explanation above 95%, confirmed the ability to represent the Factor 
(1) epidemiological indicators, pointing out, that 08 (eight) variables influenced the set of 
variables with explanation between 99.66% - 99.99%. These variables meet the group of 
access to elementary education, income, non-health public goods and services. 

It can be said that the statistical attribute commonality directly influences the result 
of the research, since it was highly present in the first Scenario (1) with 69 variables 
and in the second Scenario (2) with 16 variables. Among the variables that presented 
commonality for this study, the following can be considered as the most important: number 
of cancer preventive exams, pregnancy complication and number of deaths in the period 
corresponding to 4 years of studies. 

As for the agreed indicators, 07 (seven) are part of the indicators of service coverage 
and 26 of the epidemiological indicators. However, one can attribute a criticism to the non- 
agreement of socioeconomic indicators, which demonstrates low power of articulation in the 
field of intersectoral policy between the segments of education, public safety, environment 
and others that infers scope in the area of public health. 
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RESUMEN: El propósito de la investigación 
consistió en llevar a cabo un análisis 
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comparativo utilizando datos satelitales e 
información de estaciones meteorológicas 
para estudiar la radiación solar global en la 
ciudad de Puno. Para este fin, se tomaron 
en cuenta los datos proporcionados por el 
National Solar Radiation Data Base, las 
estaciones meteorológicas del Servicio 
Nacional de Meteorologia e Hidrologia 
del Perú en la ciudad de Puno y la 
estación de la Universidad Nacional del 
Altiplano. El análisis se centró en realizar 
comparaciones y establecer relaciones de 
corrección basadas en los datos medidos 
en la ciudad de Puno. Se consideraron tanto 
correlaciones lineales como coeficientes 
de correlación para estimar los datos de 
radiación solar global para la ciudad de 
Puno. Los resultados de la investigación 
indican que es viable obtener datos de 
radiación solar global confiables a partir de 
las bases de datos disponibles en la web, 
así como mediante el uso de estaciones 
meteorológicas. Además, se logró 
establecer relaciones lineales de corrección 
al considerar las mediciones de radiación 
solar global efectuadas por las estaciones 
meteorológicas. 

PALABRAS CLAVE: Radiación Solar, 
ciudad de Puno, correlación, estación 
meteorológica. 
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ABSTRACT: The purpose of the research was to carry out a comparative analysis using 
satellite data and information from weather stations to study global solar radiation in the city 
of Puno. For this purpose, the data provided by the National Solar Radiation Data Base, the 
meteorological stations of the National Service of Meteorology and Hydrology of Peru in the 
city of Puno and the station of the National University of the Altiplano were taken into account. 
The analysis focused on making comparisons and establishing correction relationships 
based on the data measured in the city of Puno. Both linear correlations and correlation 
coefficients were considered to estimate the global solar radiation data for the city of Puno. 
The research results indicate that it is feasible to obtain reliable global solar radiation data 
from the databases available on the web, as well as through the use of weather stations. 
In addition, it was possible to establish linear correction relationships when considering the 
measurements of global solar radiation made by the weather stations. 

KEYWORDS: Solar radiation, city of Puno, correlation, weather station. 


11 INTRODUCCIÓN 


La energia solar es la principal fuente de energia en la Tierra y la atmósfera, 
representando más del 99,97% de toda la energia sin explotar. Esta energia calienta la 
superficie terrestre y el aire, generando vientos, corrientes, evaporación, formación de 
nubes y precipitaciones. Es un impulsor fundamental de diversos procesos atmosféricos y 
terrestres, desempefiando un papel crucial en la dinámica del clima y sustentando la vida en 
nuestro planeta. (Blal et al., 2020). Cuando la radiación solar global atraviesa la atmósfera 
terrestre, experimenta interacciones complejas con diversas partículas y moléculas 
presentes en ella, como partículas de polvo, moléculas de gas, ozono y vapor de agua. 
Estas interacciones incluyen fenómenos como la reflexión, la disipación y la absorción de la 
radiación solar por parte de estas partículas y moléculas. (Cervantes-Osornio et al., 2017). 
La cantidad de atenuación de la radiación solar en un momento y lugar determinado está 
determinada por la composición de la atmósfera y la longitud del trayecto al que se transmite 
la radiación solar, es por ello la necesidad de disponer con registros de radiación solar y 
estos cobra importancia principalmente debido al incremento en aplicaciones de la energia 
solar (Sindhu et al., 2016), en áreas de la ingeniería así como en: modelación del clima y 
tiempo, monitoreo de crecimiento en los cultivos y control de enfermedades. Obtener estos 
datos confiables de radiación solar requiere mediciones sistemáticas (Mubiru, 2008). 

La radiación solar ultravioleta (UV) tiene efectos nocivos sobre la salud de los seres 
humanos. En la piel la radiación UV causando envejecimiento prematuro, reacciones de tipo 
alérgicas, como foto-queratitis (Sordo y Gutiérrez, 2013), câncer en la piel (Levels, 2017). En 
los ojos causa cataratas (Gallagher y Lee, 2006) y afecta al sistema inmunológico (Norval, 
2006). El efecto más común de la sobreexposición a la radiación solar es la quemadura 
solar o eritema, que aparece 10-16 horas después de la exposición al sol (Tanskanen et al., 
2007). Los efectos en la naturaleza como en la agricultura y la vida acuática se presenta 
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en las plantas y peces (Sweet et al., 2012). A la vegetación afectando en su crecimiento 
(Zhang et al., 2014), en los animales produce quemaduras en su piel (Pahkala et al., 2003) 
but the effects on later life-stages (metamorphosis, en la conservación y durabilidad de 
materiales como el plásticos y la pintura (Verbeek et al., 2011). 

Mayhua (2019), estudia el comportamiento de la radiación solar global para la 
ciudad de Arequipa, Perú, considerando el modelamiento matemático empírico obtenido de 
la correlación de datos de índice ultravioleta y radiación solar global obtenidos del Servicio 
Nacional de Meteorologia e Hidrologia del Perú (SENAMHI) de la estación “la Pampilla”, 
logran validar la ecuación para la estimación del índice ultravioleta. 

Antón y otros (2017) presentaron un estudio que se centra en el análisis de la radiación 
solar global diaria determinada a partir de la insolación solar en Madrid, Espafia, entre 
1887-1950. También se analiza la relación entre la radiación solar global y la nubosidad, 
así mismo la reconstrucción de la radiación solar global a largo plazo fue con tendencia 
negativa de 2 -13KJ/m por ahÃio en el periodo 1887 a 1915 y con tendencia positiva de 2 
+13KJ/m por afio en el periodo 1916 a 1950, ambos estadísticamente significativos en el 
nivel de confianza del 95%. 

Ayodele y Ogunjuyigbe (2015) proponen la distribución de probabilidad del índice 
de claridad para la predicción de la radiación solar global. El índice de claridad se obtiene 
a partir de la radiación solar global de afios pasados. Para validar el método propuesto, 
se analizaron datos durante ocho afios (2000-2007) en Ibadan, Nigeria. Los resultados 
mostraron que la distribución logística proporciona el mejor ajuste para el índice de claridad 
y el método. 

Por ello nos propusimos como objetivo realizar un análisis de la radiación solar 
global en la ciudad de Puno. 


21 MATERIALES Y MÉTODOS 


2.1 Datos de radiación solar en la ciudad de Puno 


Se utiliza datos de radiación solar medidos por la estación meteorológica DAVIS 
Vantage PRO 2.0 Plus mediante el sensor radiación solar (Piranómetro), el cual se muestra 
en lafigura 1, en intervalos de 1hora en unidades de w/m?2. La estación meteorológica está 
instalada en una zona libre de sombras y obstáculos en el pabellón de ingeniería civil de 
la Universidad Nacional del Altiplano - Puno (UNAP) en las coordenadas latitud -15,4890, 
longitud -70,1515 y una altitud de 3,825msnm. 
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(b) 


Figura 1. (a) Estación meteorológica DAVIS vantage PRO 2.0 Plus (b) Piranómetro 


También se considera la estación meteorológica de SENAMHI ubicada en la 
coordenadas latitud -15.826250 longitud -70.012083, el cual se ubica a 747.90m de la 
estación meteorológica de UNAP como se muestra en la figura 2, el cual se ubica a orillas 


del lago Titicaca. 


Figura 2. Ubicación de la estación meteorológica UNAP y SENAMHI ciudad de Puno 


2.2 Datos satelitales 

La National Solar Radiation Data Base (NSRDB) es un conjunto de datos de recursos 
solares públicos ampliamente utilizado que ha sido desarrollado y actualizado durante más 
de 20 afios para reflejar los avances en la medición y modelización de la radiación solar. 
La versión más reciente de la NSRDB utiliza productos satelitales de 30 minutos a una 
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resolución de 4x4 km que cubren el período 1998-2018 (Habte et al., 2017; Sengupta et 
al., 2018). 


2.3 Estadísticos de evaluación 


Error estândar promedio o raíz cuadrada del cuadrado medio del error (RMSE) 


n —ro2b/2 
RMSE = [nto Rsi) | (1) 
El error medio (MBE), 
n Ro 
MBE — eis: Rsi) (2) 
Coeficiente de determinación (R?) 
R2 fi [EL (Rmi-Rmi(Rsi- Rs] (3) 


- [E,(Rmi- Rm? | D(Rsi—RsiD)?] 


El coeficiente de correlación (r) 


r=vR? (4) 


Donde: R,, es el dato medido por la estación meteorológica, R,, es el dato obtenido 
del satélite, R,, es el promedio de los datos medidos por la estación meteorológica, R, es 
el promedio de los datos obtenidos del satélite y n es el número total de observaciones 
(Alexandris et al., 2006; Tabari, 2010) 


31 ANÁLISIS DE RESULTADOS 


3.1 Radiación solar del periodo 1997-2020 


La radiación solar lo obtenemos de los datos de la NSRDB (NASA), como se muestra 
en la figura 3. Sea obtenidos valores de radiación solar global para un plano horizontal para 
el periodo 1998 al 2020 de la NSRDB considerando las coordenadas de coordenadas latitud 
-15,4890, longitud -70,1515 para la ciudad de Puno con valores máximo de 9.806kwh/m"2/ 
día, mínimo de 3.004kwh/m"2/día y un promedio de 6.857kwh/m/2/dia, el cual se muestra 
enla Figura 3, se observa un comportamiento senoidal similar a los obtenidos por (Jimenez 
et al., 2016) para la ciudad de Turcuman, Argentina, y (Antón et al., 2017) obtiene valores 
de radiación solar global para la Península Ibérica en el periodo 1986-2015, encontrando 


comportamiento estacional. 
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Figura 3. Radiación solar NASA para la ciudad de Puno, afio 1998-2020. 


3.2 Radiación solar para ciudad de Puno 


Ahora se analiza datos de SENAMHI y NASA para la ciudad de Puno, para un 
periodo de 2013- 2020 de la radiación solar global diaria en unidades kwh/m"2/dia el cual 
se muestra en la figura 5. En la figura 4 se observa datos de la radiación solar para el afio 
2020, en el cual se observa que existe una diferencia entre los valores medidos, estos 
está de acuerdo al trabajo de Nolasco et al., (2018), ahora para realizar la corrección en 
función a los datos medios por la estación de SENAMHI se realiza la correlación de los 
datos SENAMHI y NASA el cual se observa en la figura 5, datos para el periodo 2013 -2020. 
Considerando la totalidad de días nublados y parcialmente nublados y despejados, donde 
sea obtenido valores mínimos, máximo y promedio los cuales se muestran en la tabla 1. E 
unidades de kwh/m"2/dia. 


—— SENAMHI 
105 —— NASA 


Elo o] 
ad dito ad 
Ea rr e ; 


Radiacion solar (kwh/m"2/dia) 


T T T T T T 
o 50 100 150 200 250 300 350 


Tiempo (Dias) 


Figura 4. Radiación solar SENAMHI y NASA para la ciudad de Puno, afio 2020. 
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Figura. 5 Comparación de datos NASA y SENAMHI 2013 a 2020. 


Valores NASA SENAMHI 
Mínimo 2.450 1.619 
Máximo 9.170 9.194 
Promedio 6.088 6.423 


Tabla 1. Datos de valores de mínimo, máximo y promedio para la el 2020 


Al realizar la comparación de datos se obtiene una ecuación lineal de relación entre 
los datos de la NASA y SENAMHI el cual es: 


Rç= 0.9826+Ry + 0.0951 (5) 


Donde R, es la radiación solar medido por SENAMHI y R, es la radiación solar 
estimado por la NASA, esta ecuación nos va permitir estimar datos para la ciudad de Puno. 

Ahora analicemos el periodo de días del 01 de abril al 23 de julio del 2020, en el 
cual se presenta un mínimo número de días con nubosidad, sin presencia de Iluvias como 
se observa en la figura 6a, donde se observa una sobre estimación de la radiación solar 
global por los datos obtenidos de NSRDB comparados con la estación, estos resultados 
concuerdan con los trabajos de (Antón et al., 2010; Buchard et al., 2008; Nolasco et al., 
2018) donde se indican que los datos obtenidos de satélites sobre estiman los valores de 
radiación solar entre 17 a 40%, esto debido a las características climáticas del lugar (Blal 
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et al., 2020). Por ello realizamos el análisis de dispersión de datos mediante el coeficiente 


de correlación lineal encontrando r= 0.952, lo que indica del total de datos para cielos sin 


nubosidad y Iluvias el 95.2% tienen relación lineal como se muestra en la figura 6b. 


Radiación solar (kwh/m"2/dia) 


SENAMHI 


75 


7.0 


6.5 


6.0 


5.5 


5.0 


45 


4.0 


3.5 


3.0 


—— SENAMHI 
—*— NASA 


22/05/2020 


29/05/2020 — 
5/06/2020 iii 


12/06/2020 ii 


19/06/2020 
26/06/2020 


3/07/2020 — 


17/07/2020 — 


31/07/2020 


7/08/2020 — 


14/08/2020 — 


21/08/2020 — 


28/08/2020 — 


4/09/2020 — 


O 24/07/2020 A 


— 


3 10/7/2020 


2 
fo) 

Es 
O 
5) 


-— 
» 
— 


m Datos 
——— Ajuste Lineal 


O um 
Rs=0.846*Rn+1.029 E 


r=0.925 n 


(b) 


Figura 6 (a) Radiación solar global (b) Dispersión de datos NSRDB y estación 


Para realizar la corrección de datos obtenidos satelitales por medio de la web según 


los autores (Jacobs et al., 2004; Mazorra Aguiar et al., 2019; Nolasco et al., 2018) se 
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considera la ecuación de la regresión lineal de los datos NSRDB y estación para la ciudad 
de Puno, se obtiene la ecuación 


Rs = 0.846 + Ra + 1.029 (6) 


Donde Rs, es la radiación solar global ajustado para la ciudad de Puno, Ra es la 
radiación solar global obtenido de NASA ambos en unidades de kwh/m"2/dia. 

Promedio de la radiación solar global para la ciudad de Puno por meses y horas para 
el periodo 2013 a 2020 de los datos medidos por la estación meteorológica de SENAMHI, se 
muestra en la figura 7, considerando la totalidad de días nublados, parcialmente nublados 
y despejados. 


Figura 7. Promedios horarios de la radiación solar global en la ciudad de Puno 2020. 


En la figura 8 se observa radiación solar global promedio mensual del periodo 
2010-2018 en el cual se observa que los avalores mínimos se presentan en la estación 
de invierno y los valores máximo se presenta en la estación de verano lo cual está de 
acuerdo a la distancia del sol y la tierra según (lgbal, 1983) y el comportamiento de los 
valores promedio mensuales es estacional y cíclico de los valores de radiación solar global 
promedio mensual. 
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Figura 8. Radiación solar global promedio mensual del periodo 2010-2018 


41 CONCLUSIONES 


Se ha llegado a la conclusión de que es factible obtener datos confiables sobre 
radiación solar a partir de las bases de datos proporcionados por el National Solar Radiation 
Data Base y las estaciones meteorológicas UNAP y SENAMHI, y estableciendo relaciones 
lineales de corrección basadas en mediciones de estaciones meteorológicas que miden 
la radiación solar global. Sin embargo, una limitación importante para estimar la radiación 
solar global en la ciudad de Puno radica en la falta de estaciones meteorológicas que 
abarquen toda la provincia de Puno. Además, las estimaciones de radiación solar global 
proporcionadas por las bases de datos de la National Solar Radiation Data Base son para 
áreas extensas, lo que puede diferir de las mediciones puntuales realizadas por estaciones 
meteorológicas. Esta discrepancia en la escala de medición puede afectar la precisión de 
las estimaciones de radiación solar global en la ciudad de Puno. 
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RESUMO: A energia solar é uma 
modalidade de energia limpa que pode 
ser utilizada como fonte de energia 
térmica em secadores convectivos. Seu 
aproveitamento pode ser realizado através 
do uso de coletores solares, responsáveis 
por transformar a energia da radiação solar 
em energia térmica, visando o aquecimento 
do ar ambiente. O objetivo deste trabalho 
foi avaliar o desempenho térmico de um 
coletor armazenador solar. O coletor 
construído foi do tipo placa plana horizontal 
e para maximizar a absorção da radiação 
solar foram utilizados, separadamente 
uma, chapa metálica de zinco e uma 
palha de aço comercial. Medidas das 
temperaturas de bulbo seco e bulbo úmido 
e das velocidades do ar na entrada e saída 
do coletor, em intervalos de uma hora no 
período de 9 às 17 horas, foram realizadas. 
Para a avaliação da eficiência instantânea 
do equipamento, foram determinados 
inicialmente o ganho de energia útil, a 
área útil do coletor e a irradiação global. 
Os resultados demonstraram que a maior 
temperatura obtida na saída coletor foi 48 
ºC, utilizando a configuração com a palha 
de aço. A irradiação global média obtida foi 
na faixa de 124 a 636 W/m?. Foi observado 
que não houve uma grande variação entre 
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as velocidades do ar na entrada e saída do coletor, como também entre as vazões mássicas. 
A configuração com a chapa metálica de zinco apresentou uma eficiência média de 52 %, 
pouco menor que a com palha de aço que foi 53 %, não havendo diferença significativa (p < 
0,05) entre elas. Assim, a utilização da palha de aço como material absorvente se demonstrou 
mais aplicável do que a chapa metálica de zinco, por ser mais acessível e possuir menor 
custo de aquisição. 

PALAVRAS-CHAVE: Coletor armazenador solar horizontal, secagem solar, eficiência 
térmica. 


THERMAL PERFORMANCE OF A HORIZONTAL SOLAR STORAGE 
COLLECTOR 

ABSTRACT: Solar energy is a clean energy that can be used as a source of thermal energy 
in convective dryers. lts use can be carried out through solar collectors, converting solar 
energy into thermal energy, aiming at heating the ambient air. The objective of this work was 
to evaluate the thermal performance of a solar storage collector. The collector built was a 
horizontal flat plate type and to maximize the absorption of solar radiation, a zinc metal sheet 
and a commercial steel wool were used separately. Measurements of dry bulb and wet bulb 
temperatures and air velocities at the collector inlet and outlet, at one-hour intervals from 9 
am to 5 pm, were performed. To evaluate the instantaneous efficiency of the equipment, the 
useful energy gain, the useful area of the collector and the global irradiation were initially 
determined. The results showed that the highest temperature obtained at the collector outlet 
was 48 ºC, using the steel wool configuration. The average global irradiation obtained was in 
the range of 124 to 636 W/m?. It was observed that there was not a great variation between 
the air velocities at the collector inlet and outlet, as well as between the mass flows. The 
configuration with zinc metal sheet showed an average efficiency of 52%, slightly lower than 
that with steel wool, which was 53%, with no significant difference (p < 0.05) between them. 
Thus, the use of steel wool as an absorbent material proved to be more applicable than zinc 
metal sheet, as it is more accessible and has a lower acquisition cost. 

KEYWORDS: Horizontal solar collector, solar drying, thermal efficiency. 


11 INTRODUÇÃO 


A secagem é uma alternativa amplamente utilizada para a redução de perdas e 
preservação da qualidade de produtos agrícolas, pois remove a sua umidade através da 
transferência de calor e massa (OLIVEIRA, 2018). Na operação de secagem, utilizam-se, 
normalmente, secadores convencionais (convectivos), tendo como fonte de aquecimento 
do ar o uso de resistências elétricas e de outros dispositivos resultantes da queima de 
diversos tipos de combustíveis, tais como: lenha, resíduos agrícolas e combustíveis fósseis 
(gás) (PRECCl et al., 2001). Esses dispositivos exigem elevado investimento e altos custos 
de funcionamento e manutenção, além de trazerem problemas ao meio ambiente. Como a 
radiação solar no Brasil é abundante, em quase todo o ano, a secagem solar se torna uma 
alternativa adequada, reduzindo custos de processamento e possíveis impactos ambientais 
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(SILVEIRA, 2016). 

A energia solar é uma modalidade de energia limpa que pode ser utilizada como 
fonte de energia térmica, para aquecimento de ambientes e de fluidos (SILVA, 2015). No 
Brasil, esta possui alto potencial de exploração, pois a irradiação solar global incidente 
em qualquer região do território brasileiro varia de 4.200 a 6.700 kWh/m?/ano, superior às 
verificadas em outros países que utilizam bastante a energia solar (PEREIRA et al., 2006). 

A Região Nordeste possui o maior potencial de energia solar disponível no Brasil, 
com uma radiação global média estimada em 5,9 kWh/m?, de acordo com o Atlas Brasileiro 
de Energia Solar (SILVA, 2015). Essa região também apresenta elevados valores anuais 
de temperatura do ar (20 a 28 “C), decorrentes, exatamente, da sua alta incidência de 
radiação solar (PEREIRA et al., 2006). 

Na região semiárida do estado da Paraíba, os pequenos produtores rurais 
sofrem com a escassez dos recursos, e por esse motivo se afiliam às cooperativas, que 
compram antecipadamente a sua produção para posterior distribuição nos mercados da 
circunvizinhança (SCHMIDT, 2008). Para que a agricultura familiar não seja prejudicada 
pelas altas temperaturas da região, seus produtores devem encontrar novos meios para 
comercializar seus produtos e de melhorar sua renda. Uma possível alternativa é a secagem 
de uma parte da produção, através da radiação solar, pois o semiárido paraibano possui 
elevada incidência solar, quando comparada a outros estados. 

A radiação solar pode ser utilizada diretamente como fonte de energia térmica 
para aquecimento de fluidos. Seu aproveitamento pode ser realizado através do uso de 
coletores ou concentradores solares (COSTA, 2010; FALCÃO FILHO, 2011). O coletor solar 
pode ser definido como um tipo de trocador de calor, responsável por transformar a energia 
da radiação solar em energia térmica, visando o aquecimento do ar ambiente (CARVALHO 
& SIONEK, 2015). Um modelo simples de coletor é aquele constituído por uma câmara 
plana ou suporte (geralmente de madeira), um tampo de vidro, intercalados por uma 
chapa metálica (coletora). A radiação atravessa o vidro, incide na placa coletora e passa a 
aquecer o ar que atravessa o sistema, no qual ocorre uma diferença de temperatura e de 
taxa mássica entre as suas correntes de entrada e saída do coletor, resultando em uma 
descarga de ar aquecido desse equipamento (BECKER, 2014). 

O presente trabalho teve como objetivo principal a construção de um coletor 
armazenador solar plano horizontal de custo mais acessível a pequenos produtores rurais 
das áreas mais secas do Brasil, que poderão realizar a secagem de sua produção agrícola 
facilitando a comercialização de seus produtos e aumentando sua renda familiar. 
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21 MATERIAL E MÉTODOS 


Realização dos experimentos 


Os experimentos foram conduzidos no Laboratório de Engenharia de Alimentos 
(LEA) do Departamento de Engenharia de Alimentos (DEA/CT/UFPB), João Pessoa - PB. 
O clima da região onde foram realizados os ensaios apresenta temperatura média anual de 
25 ºC e precipitação média anual de 1.900 mm (INMET, 2018). 


Desenvolvimento do projeto do coletor armazenador solar 


O coletor solar construído foi do tipo placa plana horizontal. Esta escolha se deve ao 
fato de se optar por uma convecção forçada do ar. Estudos comprovam que a velocidade do 
ar possui uma maior influência no início do processo de secagem sobre a menor temperatura, 
não exercendo influência significativa na maior temperatura de secagem (ARAÚJO, 2005; 
EL AOUAR, 2002; AZOUBEL, et al., 2009), o que justifica a utilização de uma convecção 
forçada do ar. O referido coletor foi constituído de uma caixa em aglomerado de madeira, 
com dimensões de 1,00 x 0,50 x 0,15 m; 4 separadores internos (290,0 x 10,0 x 150,0 
mm), também em aglomerado de madeira com a finalidade de aumentar o percurso do ar 
para aproveitar melhor a área de transferência de calor; um soprador para a movimentação 
forçada do ar; e um tampo de vidro. 

Toda a estrutura foi fixada por meio de parafusos e devidamente vedada com silicone. 
Este dimensionamento foi determinado de acordo com recomendações de Bezerra (2001), 
onde para cada metro cúbico da câmara de secagem é necessária uma área de 8 m? de 
coletor solar, formando uma proporção de 1 m3:8 m2. 


O desenho esquemático do coletor solar está apresentado nas Figuras 1 e 2 a seguir. 


Aa 
“ 


Ê 
Ê 
| 


| Coletor com chapa metálica de zinco Coletor com palha de aço comercial 


Figura 1 — Vista superior do coletor solar com (a) chapa metálica de zinco e com (b) palha de aço 
comercial. 


Fonte: Autor (2023). 
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100cm 


Figura 2 — Vista lateral do coletor solar. 


Fonte: Autor (2023). 


Para maximizar a absorção da radiação solar foram utilizados separadamente dois 
elementos absorvedores distintos: uma chapa metálica de zinco (CM) e uma palha de aço 
comercial (PA). Esses elementos foram cobertos por tinta de cor preta, para promover uma 
maior absorção. As dimensões dessas duas últimas partes constituintes do coletor foram 
acertadas no momento de sua construção. A espessura de PA utilizada foi de 2,5 cm. A 
convecção forçada de ar foi obtida a partir de um ventilador axial de baixo custo, conforme 
ilustrado nas Figuras 1 e 2. 


Análise térmica do coletor 


Foram realizados testes na área externa do Laboratório de Engenharia de Alimentos 
(LEA), com o intuito de obter dados suficientes para o cálculo da eficiência instantânea, 
bem como caracterizar as condições psicrométricas do ar ambiente e do ar que percorreu 
o coletor. Foi utilizado um termômetro de mercúrio, para obter temperaturas de bulbo seco 
e úmido, e um anemômetro de fio quente da marca AKSO modelo AK833 para verificar as 
velocidades na entrada e saída do coletor. 

Em todos os testes, o coletor ficou na posição horizontal e os dados de temperatura 
e radiação solar foram registrados com intervalos de uma hora. Os valores de temperaturas 
de bulbo seco e úmido coletadas na entrada e saída do coletor foram suficientes para obter 
as propriedades termodinâmicas do ar, com o auxílio do Sofware Grapsif desenvolvido por 
Melo et al. (2004). 


Desempenho do coletor solar 


Após a coleta dos dados experimentais, foi realizada a avaliação do desempenho 
do coletor solar, seguindo a metodologia proposta por Silveira (2016). O desempenho do 
coletor foi determinado pela eficiência (1), obtida pela equação 1 (DUFFIE & BECKMAN, 
2013). 
ni= e eq.1 


= AcçX Gr 


O cálculo da eficiência foi realizado para ambas as configurações, contudo, para 
melhor comparação entre elas optou-se por usar uma média das eficiências obtidas nos 
ensaios com a chapa metálica de zinco para comparar com as eficiências adquiridas no 


ensaio com a palha de aço. Para poder calcular a eficiência, precisou-se saber inicialmente 
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o ganho de energia útil (Qu) (vide equação 2) (KALOGIROU, 2014), a área útil do coletor 
(Aq), encontrada a partir da medição do comprimento e largura úteis do coletor, por onde o 
ar circula, e também o valor da irradiação global (1G). 


Qu =mxcCpX(Tes,— Tas) eq.2 


Para a obtenção dos dados da irradiação global nos horários da medição (de hora 
em hora), foi necessário um levantamento no site do Instituto Nacional de Meteorologia 
(INMET). O site disponibiliza dados de temperatura máxima e mínima, ponto de orvalho, 
umidade relativa, como também a radiação e a quantidade de chuva em um determinado 
local a cada hora. Como o valor disponível era de radiação (R) em unidade de energia (kJ/ 
m?), e para o cálculo da eficiência instantânea é necessária a irradiação (IG) em unidade 
de potência (W/m?), foi feita a conversão utilizando a equação 8. 


R x1000 
3600 


IG = 


eq.3 


Multiplicou-se por mil para transformar de kJ para J, e dividiu-se por 3600 para 
converter horas em segundos, e passar a ser J/s que é o mesmo que W. Foi feita essa 
conversão para cada valor coletado. No site o horário era dado em UTC, por isso foi preciso 
fazer um ajuste para o estado da Paraíba, adicionando três horas ao horário informado no 
site. 

O ganho de energia útil foi determinado a partir da equação 2, que relaciona as 
temperaturas de entrada (TBSe) e saída (TBSs) do coletor, o calor específico do ar (cp), 
obtido segundo Singh (2009), e a vazão mássica de ar (Mm), calculada pela equação 4. 
Essa última grandeza, relaciona a média entre as velocidades de saída e entrada (v,), a 
área da seção transversal (Aq), obtida a partir do diâmetro da tubulação onde é coletada a 
velocidade de saída. 


mM = AçX PaX Va eq.4 


A massa específica do ar (p,) foi encontrada pela equação 5, desenvolvida por 
Monteiro et al. (2010), que depende da temperatura média de entrada e saída do coletor 
(To). 


m 


— 3534 
Pa = Tm+273,15 


eq.5 


Tratamento estatístico dos dados experimentais 


Os dados obtidos experimentalmente foram analisados por meio do pacote estatístico 
STATISTICA 7.0, mediante Análise de Variância (ANOVA) sobre a eficiência instantânea do 


coletor, para ambas as configurações mencionadas anteriormente. 
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31 RESULTADOS E DISCUSSÃO 


Desenvolvimento do projeto do coletor 


O coletor construído foi baseado na metodologia desenvolvida por Silva (2013), 
que tinha como intuito construir um coletor solar de baixo custo, para aplica-lo na secagem 
de produtos alimentícios. Outra metodologia que inspirou esse trabalho foi a de Oliveira 
(2018), que desenvolveu um coletor solar acoplado a um secador híbrido. 


Na Figura 3 é possível ver ambas as configurações testadas. 


Figura 3 — Coletor solar com as configurações (a) CM e (b) PA. 
Fonte: Autor (2023). 


Análise térmica do coletor 


A análise térmica do coletor foi realizada a partir das medidas de temperatura de 
bulbo seco na entrada e saída do mesmo, assumindo duas configurações distintas: o 
coletor com a chapa metálica de zinco (CM) e com a palha de aço comercial (PA). 

A Tabela 1 apresenta os dados de temperaturas na entrada e saída do coletor e 
irradiação global, para ambas as condições mencionadas anteriormente, ao longo do dia. 


CM PA 
Ei cê) (6) (im cê) (6 (Wim 

9 33 a” 458 30 aa 487 
10 34 39 606 aa 36 463 
11 34 39 639 34 39 589 
12 33 38 591 38 48 595 
13 35 39 434 35 41 838 
14 35 40 555 32 35 432 
15 34 36 455 33 35 335 


Estudos multidisciplinares: Ciências exatas e da terra e engenharias Capítulo 3 


43 


16 30 31 301 32 35 284 
17 30 30 98 31 31 150 


Tabela 1 — Análise térmica do coletor solar, comparando as duas configurações utilizadas (CM e PA) 
em relação às temperaturas de bulbo seco de entrada e saída e irradiação global. 


Fonte: Autor (20283). 


Observando os dados da Tabela 1, percebe-se que, na primeira hora do 
experimento, já houve uma variação de 4ºC na temperatura da configuração CM, enquanto 
na PA aumentou 3 “C na temperatura de entrada do coletor. No horário de 12 horas, a 
temperatura de bulbo seco na saída da configuração CM teve uma queda com relação a 
medição anterior, ficando 38 ºC, diferente da configuração PA, que já iniciou a medição 
com essa temperatura na entrada, chegando a alcançar, na saída, a maior temperatura de 
todas as medições (48 “C). No fim da tarde, as temperaturas das duas configurações foram 
iguais na entrada e na saída, pois a posição do coletor nesse momento era na sombra. 
Além disso, foi observado que coletor apresentou características de armazenador para 
a configuração PA, quando às 16 horas, mantendo praticamente a mesma variação de 
temperatura em relação às 15h. 

A irradiação global média variou de 124 a 636 W/m?, aproximadamente, em uma 
medição que se iniciou de 9 horas e foi até às 17 horas, conforme ilustrado na Figura 4. 
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Figura 4 — Irradiação global média em função da hora do dia. 


Fonte: Autor (2023). 


A curva apresentada na Figura 4 mostra que no período de 9 a 11 horas houve um 
crescimento no valor da irradiação, reduzindo um pouco ao meio dia. O período do dia 
com maior irradiação média global foi de 13 horas, com um valor de 636,22 W/m?. Oliveira 
(2018) obteve, em junho na cidade de Seropédica-RJ, uma irradiação máxima de 593,67 
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W/m?, no período de meio dia, valor bastante semelhante ao encontrado nesse trabalho 
(593,01 W/m?). Já Silva (2013), obteve um valor médio de 611,60 W/m?, na cidade de 
Coimbra em Portugal, em um dia nublado do mês de setembro. Após alcançar o pico, o 


valor cai bruscamente, e mantém esse comportamento até o fim do dia. 


Desempenho do coletor solar 


A Figura 5 mostra a relação entre a temperatura de bulbo seco de saída e a irradiação 
solar global, comparando as duas configurações utilizadas. 

Observa-se na Figura 5 que a relação entre a variação da irradiação global e a 
variação de temperatura de bulbo seco de saída apresentou comportamento próximo à 
função linear. 

O coeficiente de determinação R2 indica o quão explicativo é o modelo linear; para a 
configuração CM, seu valor foi de 0,84 e para a configuração foi de 0,69. Portanto, o melhor 
modelo linear é o da configuração CM, que explica 84 % da variância dos dados. Valores 
condizentes com um comportamento linear da relação acima mostrada foram encontrados 
por Oliveira (2018) (0,82) e por Silveira (2016) (0,98). 
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Figura 5 — Temperatura de bulbo seco de saída (ºC) para cada configuração por irradiação global 
média (W/m?). 


Fonte: Autor (2023). 


Ainda pela Figura 5, percebe-se que as equações das retas para ambas as 
configurações são, praticamente, iguais, o que confirma o observado na prática. 

De acordo com a metodologia proposta por Silveira (2016), foi possível obter os 
valores médios de velocidade, massa específica e vazão mássica do ar, os quais foram, 
respectivamente, de 3,67 m/s, 1,15 kg/mº e 0,0331 kg/s. A área da seção transversal foi de 
0,00785 m?, obtida a partir do diâmetro da tubulação de saída (0,1 m) (MONTEIRO et al., 
2010). 

A partir dos valores mencionados, foram calculadas as eficiências do coletor, para 
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ambas as configurações, a cada hora do dia, conforme a Figura 6. 
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Figura 6 — Eficiência instantânea de cada configuração por hora do dia. 


Fonte: Autor (2023). 


Verifica-se, na Figura 6, que nos primeiros horários de medição, o coletor com a 
configuração CM se mostrou mais eficiente que o com a configuração PA. Já das 11 até às 
12 horas o mais eficiente foi o coletor com a PA. A menor eficiência instantânea foi no fim 
da tarde, graças à sombra formada no local onde estava localizado o coletor. 

Hedge et al. (2015) avaliaram a eficiência para um sistema de secagem solar de 
banana, comparando os resultados para passagem do fluxo de ar entre a placa absorvedora 
e a cobertura do coletor e para a passagem do fluxo entre a placa absorvedora e uma 
camada de isolamento. Esses autores iniciaram a medição também às 9 horas e atingiram 
a eficiência máxima de 50% às 16 horas. 

Fudholi et al. (2015) obteve 28 % de eficiência do coletor solar para a secagem 
de pimenta vermelha com fluxo de massa de ar de 0,07 kg/s. Já Lingayat et al. (2017) ao 
avaliar um secador solar indireto para a secagem de banana estimaram a eficiência média 
do coletor em 31,50 %, em condições de irradiação solar média de 724 W/m?. 

Oliveira (2018) obteve um valor de 25,9 % de eficiência média global em um coletor 
solar plano acoplado em secador híbrido para secagem de produtos agrícolas. Enquanto 
isso Potdukhe e Thombre (2008) obtiveram 34 % de eficiência em um coletor na secagem 
solar de pimentões e folhas de feno-grego. Em todos esses estudos os autores encontraram 
valores de eficiência inferiores às duas configurações apresentadas nesse trabalho. 

A configuração com a chapa metálica de zinco apresentou uma eficiência média de 
52 %, pouco menor que a com palha de aço que foi 53 %, não havendo diferença significativa 
a p < 0,05. Com isso, a utilização da palha de aço como material absorvente mostrou-se, 


aparentemente, mais eficiente que a chapa metálica de zinco, além de ser mais vantajosa 
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com relação aos custos, pois uma folha de zinco de dimensões 100 x 100 x 0,5 mm custa 
em média R$ 10,18 no site do AliExpress, e para cobrir a área do coletor são necessárias 
10 folhas, totalizando um custo de R$ 101,80. Já a palha de aço pode ser encontrada em 
qualquer supermercado com o preço variando de R$1,50 a R$3,00, sendo Bombril a marca 
mais conhecida. Para a utilização no coletor, seriam necessários 4 pacotes do produto, o 
que não ultrapassaria R$12,00, tendo uma economia de aproximadamente 88 %. 


41 CONCLUSÕES 


O coletor construído apresentou algumas perdas resultantes do tipo de ventilador 
escolhido e da estrutura onde este se encontrava acoplado. A maior temperatura de bulbo 
seco na saída do coletor e a maior irradiação média global foram obtidas utilizando a 
configuração com a palha de aço, no período de 12 às 13 horas, mesmo em dias nublados. 

Foi observado durante os experimentos que não houve uma grande variação tanto 
entre as velocidades de entrada e saída do coletor quanto entre as vazões mássicas do 
ar. A configuração com a chapa metálica de zinco apresentou uma eficiência média pouco 
menor que a com palha de aço. A análise estatística mostrou que não houve diferença 
significativa entre as duas configurações estudadas, à 5% de significância. Portanto a palha 
de aço comercial pode ser utilizada como meio absorvedor ao invés da chapa metálica 
de zinco, pois é bem mais acessível tanto do ponto de vista de disponibilidade quanto de 
preço. 
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RESUMO: Simulações termoenergéticas 


foram realizadas em uma edificação 
localizada em Belém-PA | utilizando 
os softwares SketchUp, OpenStudio, 


EnergyPlus e DView para avaliar o 
comportamento térmico do espaço e 
determinar a representatividade (%) em 
que o ambiente se mantém dentro dos 
parâmetros estabelecidos pela Norma 
ISO 9241. Cinco arquivos climáticos 
diferentes foram utilizados para comparar 
o desempenho térmico do ambiente: 
TRY, IWEC, SWERA, INMET e TMY. Com 
base na análise da curva de duração, 
constatou-se que o arquivo climático TRY 
teve o melhor desempenho, pois a zona 
termal permaneceu dentro dos parâmetros 
estabelecidos pela norma ISO 9241 por 
2.960 horas, correspondendo a 33,8% do 
período anual. 
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climáticos. Comportamento térmico. 


STUDY OF DURATION CURVES FOR THERMAL AND STRUCTURAL 
BEHAVIOR EVALUATION USING DIFFERENT CLIMATE FILES FROM THE CITY 
OF BELÉM-PA 


ABSTRACT: Thermoenergetic simulations were carried out in a building located in Belém- 
PA using SketchUp, OpenStudio, EnergyPlus and DView software to evaluate the thermal 
behavior of the space and determine the representativeness (%) in which the environment 
remains within the parameters established by ISO 9241 Standard Five different climate files 
were used to compare the thermal performance of the environment: TRY, IWEC, SWERA, 
INMET and TMY. Based on the analysis of the duration curve, it was found that the TRY 
climate file had the best performance, as the thermal zone remained within the parameters 
established by the ISO 9241 standard for 2.960 hours, corresponding to 33.8% of the annual 
period. 

KEYWORDS: Simulations. Climate files. Thermal behavior. 


11 INTRODUÇÃO 


O comportamento térmico de um edifício é uma questão crítica para a eficiência 
energética. É importante entender como o clima e outros fatores externos afetam o 
desempenho termoenergético de um edifício para que se possa tomar decisões informadas 
sobre como melhorar sua eficiência. 

A análise do comportamento térmico em um edifício pode ser feita por meio de 
modelagem de simulação numérica, que permite a previsão do desempenho em diferentes 
condições climáticas. 

Um dos métodos para avaliar o comportamento térmico é através do uso do gráfico 
de curva de duração. Este gráfico fornece uma visualização de como a temperatura 
interna do ar varia ao longo do tempo, permitindo uma análise aprofundada das condições 
climáticas e da eficiência termoenergética. 

O objetivo deste artigo é avaliar o comportamento térmico do ar em uma edificação 
por meio da análise da curva de duração utilizando os softwares SketchUp, OpenStudio, 
EnergyPlus e DView. Além disso, serão avaliados os impactos de diferentes arquivos 
climáticos na análise, os resultados são comparados com os parâmetros da norma 
reguladora ISO 9241 (2011). 

Ao compreender o comportamento térmico do ar de um edifício por meio da curva 
de duração, pode-se tomar decisões precisas sobre como gerenciar a energia de maneira 
mais eficiente e sustentável. Ao utilizar os softwares de simulação termoenergética pode- 
se ter uma visão clara e precisa sobre como melhorar a eficiência termoenergética de um 
ambiente. 

Ao utilizar arquivos climáticos diferentes, pode-se analisar como as condições 
climáticas afetam o comportamento térmico do ar em um edifício. e.g. se uma edificação for 
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simulada com dados climáticos de uma região com clima quente e úmido, pode-se prever 
como o comportamento térmico deste ambiente seria diferente em uma região com clima 
frio e seco. Isso permite avaliar as opções de gerenciamento de energia mais eficientes para 
diferentes condições climáticas e determinar qual é a melhor solução para cada situação. 
Em resumo, a análise do comportamento térmico de um ambiente estrutural utilizando 
arquivos climáticos diferentes via softwares de simulação termoenergética se torna uma 
ferramenta necessária para avaliar a eficiência térmica de um ambiente, baseando-se em 


normas regulamentadoras. 


21 REVISÃO BIBLIOGRÁFICA 


De acordo com Frota e Schiffer (2003), o conforto térmico é a sensação que 
o corpo humano experimenta ao descarregar calor produzido pelo metabolismo sem 
precisar usar mecanismos de regulação térmica. Bittencourt e Cândido (2008) afirmam 
que conforto térmico é a condição de bem-estar psicológico quanto às condições térmicas 
de um ambiente que mantém a estabilidade corporal. A norma americana ASHRAE 55 
(2004) define conforto térmico como a situação em que uma pessoa está satisfeita com 
seu entorno. Alguns autores também usam o termo conforto higrotérmico, mas a literatura 
mais recente refere-se ao conforto térmico. A ASHRAE 55 (2004) estabelece que, para se 
ter conforto térmico, pelo menos 80% dos ocupantes devem expressar satisfação com as 
condições térmicas. 

Akutsu, Sato e Pedroso (1987) enfatizam que a satisfação com o ambiente 
térmico é uma resposta que depende de diversos fatores, incluindo condições ambientais 
(temperatura, umidade, velocidade do ar e radiação) e pessoais (vestimenta e atividade). 
Romero (2013) ressalta a importância de compreender como as variáveis ambientais 
afetam a percepção térmica humana. Para garantir o equilíbrio entre o homem e seu 
meio ambiente, segundo Rivero (1986), é preciso garantir o conforto térmico em qualquer 
ambiente, seja interno ou externo. 

Diante da necessidade de definir os parâmetros de conforto térmico, foram realizados 
inúmeros estudos bioclimáticos para equacionar a satisfação humana quanto às variáveis 


climáticas. 


2.1 Norma regulamentadora ISO 9241 (2011) 

A norma recomenda que as temperaturas do ambiente sejam mantidas dentro de 
certos limites para garantir o conforto dos usuários, sendo eles: a temperatura ideal, no 
verão entre 20 “C e 24 “C e no inverno entre 23 “C e 26 ºC, com umidade relativa do ar 
entre 40% e 80%. 

2.2 Arquivos climáticos 


São conjuntos de dados que descrevem as condições climáticas de uma determinada 
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área, incluindo informações sobre temperatura, umidade, velocidade do ar, precipitação e 
outros fatores. Nesta pesquisa foram utilizados os seguintes arquivos climáticos da cidade 
de Belém-PA: International Weather for Energy Calculations (IWEC), Solar and Wind 
Energy Resource Assessment (SWERA), Instituto Nacional de Meteorologia (INMET), Test 
Reference Year (TRY) e Typical Meteorological Year (TMY). 


2.3 Softwares de simulação termoenergética 


São programas que simulam as condições climáticas e as cargas térmicas em 
edifícios e outros tipos de estruturas. Nesta pesquisa foram utilizadas as seguintes 
ferramentas computacionais: O software de modelagem 3D SketchUp (Figura 1) foi utilizado 
para modelar o ambiente, enquanto o plug-in OpenStudio (Figura 2) serviu como conexão 
entre as diferentes ferramentas de software. O EnergyPlus (Figura 3) foi utilizado para as 
simulações térmicas, e os resultados foram visualizados por meio do DView (Figura 4), que 
forneceu os gráficos que serão apresentados na seção de resultados e discussões. 
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Figura 1. Interface do software SketchUp. Fonte: Autores, (2023). 


Estudos multidisciplinares: Ciências exatas e da terra e engenharias Capítulo 4 


52 


É Untitled" 
File Preferences Components & Measures Help 


“12] Weather File & Design Days ( Utility Bits | MyMedel | Library Edit 
Weather File | Set Weather File Select Year by: 


Name: 


Elo 


Latitude: (9) Calendar Year 2022 


Longitude: O First Day of Year [UseWieatherfile 4) 
Elevation: 
Time Zone: Daylight Savings Time: (| 
Download weather files at ww energyplus.net/wenther 


Starts 


(O Define by Day of The Week And Month +) [Monday 


Measure Tags (Optional): (O Define by Date 
ASHRAE Climate Zone 


Ends 


CEC Climate Zone 


(O Define by Day of The Week And Month [First +) [Monday 


(9) Define by Date 1071/22 


Design Days | Import From DDY | 


Design Days 


gogaguddgagdat 


&; EP-Launch 
File Edit View Help 


Group of Input Files | History | Utilities | 


— Input File 


| CAUsersADavid SousatDesktop+Sem nomeidf. idf 2 | 
Browse... | Edit - Text Editor | Edit - IDF Editor | 


weather File 


ESC&T - UFPASTCCASoftwarestOS - Arquivo climático BelémiBR4, BelemVal.de.Caes.821930 SWERA.epv » 


Browse... | 


View Results 
Ê Tables | Errors DE IN | ELDMP | BND | Bsmt Out | Bsmt CSV 
Meters | RDD DE OUT | DFDMP. | DBG | Bsmt | EDD | 
WVariables MDD MAP | Screen | SLN | Bsmt Audit | Table xML | 
EIO | MTD EXPIDF | SHD | ESO | Slab Out | PerfLog.csv 
5vG Z5Z EPMIDF | vVAML | MTR | Slab | 
DxE | 552 EPMDET | Budit | Proc CSV Slab Er | 


EnergyPlus 22.2.0 


Figura 3. Interface do software EnergyPlus. Fonte: Autores, (2023). 
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Figura 4. Interface do software DView. Fonte: Autores, (2023). 


31 CARACTERIZAÇÃO DO TERRITÓRIO 


Belém-PA é a capital do estado brasileiro do Pará, situada na região Norte do país. 
A caracterização do território é uma etapa fundamental para compreender o desempenho 
térmico de ambientes. Isso ocorre porque a localização geográfica, o clima e outros fatores 
ambientais podem ter um grande impacto na temperatura, umidade e outros aspectos do 
ambiente. 

A Tabela 1 caracteriza as informações da Cidade de Belém-PA. 


POPULAÇÃO 
População estimada 1.506.420 [2021] 
População no último censo 1.393.399 [2010] 
Densidade demográfica 1.315,26 habitantes/km? [2010] 
MEIO AMBIENTE 
Área urbanizada 147,35 km? [2019] 
Esgotamento adequado 67,9 % [2010] 
Arborização de vias públicas 22,3 % [2010] 
TERRITÓRIO 
Área territorial 1.059,466 km? [2021] 


Tabela 1. Dados da cidade de Belém-PA. Fonte: IBGE — Cidades e Estados, (2023). 


A Tabela 2 caracteriza a localização geográfica da Cidade de Belém-PA. 
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Município Localização Geográfica Altitude (m) 
Latitude sul Longitude W-Gr 
Belém 01º27'20” 48º"30"15” 
Tabela 2 — Localização geográfica da cidade de Belém-PA. Fonte: CODEM, (2011). 


22 


A Figura 5 mostra o mapa da Região Metropolitana de Belém. 
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Figura 5. Região metropolitana de Belém-PA. Fonte: Google Earth, (20283). 

O Mapa de Zoneamento Bioclimático é uma ferramenta que divide o território 
brasileiro em zonas bioclimáticas, levando em consideração características climáticas e 
geográficas, como temperatura, umidade, altitude, latitude, entre outras. O objetivo do 
mapa é auxiliar na concepção de projetos de arquitetura e urbanismo mais adequados ao 
clima e condições ambientais de cada região. 

O mapa é uma importante ferramenta para orientar o desenvolvimento de projetos 
arquitetônicos e urbanísticos mais eficientes e sustentáveis, que considerem as condições 
climáticas de cada região do país. Com isso, é possível proporcionar mais conforto térmico 
e melhor qualidade de vida aos usuários dos espaços construídos, além de reduzir o 
consumo de energia e minimizar os impactos ambientais. 

O zoneamento bioclimático brasileiro compreende oito diferentes zonas, conforme 


indicado na Figura 6. 
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Figura 6. Zoneamento bioclimático brasileiro. Fonte: NBR 15220-3, (2005). 


41 MATERIAIS E MÉTODOS 


Realizou-se simulações numéricas utilizando os softwares SketchUp (v21.0.391), 
OpensStudio (v3.5.1), EnergyPlus (v22.2.0) e DView (v1.2.0) juntamente dos arquivos: 
IWEC, SWERA, INMET, TRY e TMY da cidade de Belém-PA para analisar o comportamento 
térmico do ar em um edifício por meio da análise da curva de duração, com isso, exibir 
a representatividade (%) em que o ambiente estrutural fica dentro dos parâmetros 
estabelecidos pela norma ISO 9241 (2011), avaliado a diferença de resultados com o uso 
de diferentes arquivos climáticos utilizados. Por fim, é mostrado um gráfico com a média 
dos resultados dos arquivos climáticos gerando um mapa de calor, mostrando os pontos 
quentes e frios dentro do edifício. 

A partir do software SketchUp, foi modelado um ambiente em 3D (Figura 7). Os 
materiais utilizados na construção foram aplicados a partir do OpenStudio, e a simulação 
realizada por meio do EnergyPlus. Os resultados foram apresentados em forma de gráficos 
através do DView. 
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Figura 7. Edifício básico modelado no software SketchUp. Fonte: Autores, (2023). 


As dimensões do edifício simulado são: 7 m de comprimento, 5 m de largura e 3 m 
de altura, com uma área total de 35 m?. 
Os materiais construtivos do edifício e suas propriedades térmicas são mostrados 


na Tabela 83. 

Material Espessura (m) de o p= (kg/ as 
Piso cerâmico 0,010 0,90 1600 920 
Argamassa 0,015 1,15 1600 1000 
Reboco 0,025 1,15 1800 1000 
Telha cerâmica 0,020 0,70 1000 920 

Porta (madeira) 0,030 0,15 614 2300 
Janela (vidro) 0,003 0,90 - - 

Bloco cerâmico 0,096 0,90 1764 920 


Tabela 3. Propriedades térmicas dos materiais. Fonte: ABNT NBR 15220-2, (2003). 


51 RESULTADOS E DISCUSSÕES 


Arquivo climático IWEC: o resultado apresentado na Figura 8 indica que a zona 


térmica do edifício em questão está fora dos parâmetros de conforto térmico na maior parte 
do ano, representando 74,5% do período anual. Isso significa que, durante 6.500 horas, a 
temperatura está acima de 26 “C e 30 horas abaixo de 24ºC, o que pode tornar o ambiente 
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desconfortável para os ocupantes. 

Por outro lado, a figura também mostra que a zona térmica do edifício está dentro 
dos parâmetros de conforto térmico durante 2.230 horas, o que representa 25,5% do 
período anual. Isso significa que há momentos em que a temperatura está dentro do 
intervalo desejado para o conforto dos ocupantes. 


O 500 1000 1500 2000 2500 3000 3500 4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 
Perfil Anual (h 


Figura 8. Perfil anual da curva de duração (IWEC). Fonte: Autores, (2022). 


Arquivo climático SWERA: O resultado apresentado na Figura 9 indica que a zona 


térmica do edifício em questão apresenta um grande desafio para o conforto térmico dos 
ocupantes, uma vez que a temperatura está acima dos parâmetros de conforto térmico 
durante a maior parte do ano, representando 81% do período anual, ou seja, durante 7.100 
horas. Demonstrando que a maioria dos ocupantes pode experimentar desconforto térmico, 
afetando negativamente a produtividade, o bem-estar e a saúde dos ocupantes. 

Além disso, a figura mostra que a zona térmica do edifício está dentro dos parâmetros 
de conforto térmico apenas durante 1.660 horas, representando 19% do período anual. 
Esse período é relativamente curto, o que pode ser insuficiente para garantir um ambiente 
confortável para os ocupantes. 
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Figura 9. Perfil anual da curva de duração (SWERA). Fonte: Autores, (2022). 


Arquivo climático INMET: O resultado apresentado na Figura 10 é semelhante ao 


apresentado na Figura 8, mostrando que a zona térmica do edifício está fora dos parâmetros 
de conforto térmico durante 6.600 horas, pois a temperatura está acima de 26 “Ce 30 horas 
abaixo de 24 ºC, o que representa 75,6% do período anual. 

Esses dados indicam que o ambiente pode ser desconfortável para os ocupantes e 
isso pode afetar negativamente sua produtividade, bem-estar e saúde. 

Além disso, a figura mostra que a zona térmica do edifício está dentro dos parâmetros 
de conforto térmico durante 2.130 horas, representando 24,4% do período anual. Embora 
seja relativamente baixo, esse valor é importante e deve ser valorizado, pois representa 


oportunidades para que os ocupantes se sintam confortáveis durante o ano. 
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Figura 10. Perfil anual da curva de duração (INMET). Fonte: Autores, (2022). 


Arquivo climático TRY: a Figura 11 indica que houve uma melhora significativa no 


conforto térmico da zona em comparação com os resultados dos demais arquivos climáticos 
apresentados. Durante 2.960 horas, a temperatura está dentro dos parâmetros de conforto 
térmico, representando 33,8% do período anual. Esse valor é o mais alto encontrado em 
todo o estudo, o que é um resultado positivo. 

Embora a maioria do período ainda esteja fora dos parâmetros de conforto térmico, 
com a temperatura acima de 26 “C por 5.800 horas, o fato de a zona estar dentro dos 
parâmetros de conforto térmico por quase um terço do ano é um avanço significativo. Isso 
significa que os ocupantes terão mais oportunidades para trabalhar e viver em um ambiente 
confortável. 

Além disso, é importante notar que as horas fora dos parâmetros de conforto térmico, 
representando 66,2% do período anual, ainda são consideravelmente altas, o que indica 
que ainda há espaço para melhorias no projeto e na operação do edifício. 

Portanto, é importante que medidas adicionais sejam adotadas para melhorar 
o conforto térmico da zona, como a implementação de sistemas de ventilação e ar- 
condicionado mais eficientes, uso de materiais de construção com propriedades térmicas 
isolantes e consideração da orientação do edifício e uso de elementos de sombreamento 
adequados. 
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Figura 11. Perfil anual da curva de duração (TRY). Fonte: Autores, (2022). 


Arquivo climático TMY: a Figura 12 indica que a zona térmica do edifício apresenta 


sua temperatura fora dos parâmetros de conforto térmico durante 85,1% do período anual. 
Durante 7.400 horas, a temperatura está acima de 26 “C e 60 horas abaixo de 24 ºC. 
A figura também mostra que a zona térmica do edifício está dentro dos parâmetros de 
conforto térmico durante 1.300 horas, representando 14,9% do período anual. Embora 
esse valor seja relativamente baixo, ele representa oportunidades para que os ocupantes 
se sintam confortáveis durante o ano. 
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Figura 12. Perfil anual da curva de duração (TMY). Fonte: Autores, (2022). 


A Tabela 4 exibe o resumo dos resultados das curvas de duração. 
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CURVAS DE DURAÇÃO 


Arquivo climático Dentro dos parâmetros Fora dos parâmetros 
(ISO 9241) (ISO 9241) 
IWEC 25,5% 74,5% 
SWERA 19,0% 81,0% 
INMET 24,4% 75,6% 
TRY 33,8% 66,2% 
TMY 14,9% 85,1% 


Tabela 4. Curvas de duração (resumo). Fonte: Autores, (2023). 


Através da análise das curvas de duração da temperatura na zona térmica do 
edifício, pode-se observar que o desempenho mais satisfatório em termos de conforto 
térmico foi registrado utilizando o arquivo climático TRY. Durante 2.960 horas, a temperatura 
permaneceu dentro dos parâmetros definidos pela norma ISO 9241, o que equivale a 
33,8% do período anual. Esse é um resultado positivo, indicando que a escolha do arquivo 
climático pode ter um grande impacto no conforto térmico de um edifício. 

Embora a zona ainda tenha ficado fora dos parâmetros de conforto térmico por 66,2% 
do período anual, o resultado obtido com o arquivo climático TRY é o mais alto encontrado 
em todo o estudo. Isso indica que a seleção cuidadosa do arquivo climático pode ser um 
fator crítico para garantir um ambiente confortável para os ocupantes do edifício. 

É importante ressaltar que, apesar do bom desempenho observado no arquivo 
climático TRY, ainda há espaço para melhorias na busca de soluções mais eficientes para 
garantir o conforto térmico da zona. Isso pode incluir a implementação de estratégias 
passivas, como o uso de sombreamento adequado e o emprego de materiais de construção 
com propriedades térmicas isolantes, além do uso de sistemas de ventilação e ar- 
condicionado mais eficientes. 

Por fim, a escolha do arquivo climático é um fator crítico para garantir o conforto 
térmico de um edifício, e a seleção cuidadosa do arquivo pode ter um grande impacto na 
satisfação dos ocupantes. 


61 CONCLUSÃO 


Pode-se afirmar que a análise da curva de duração é uma técnica de grande 
importância em estudos relacionados à eficiência energética e conforto térmico em 
edificações. Os resultados obtidos indicam que a utilização do arquivo climático TRY 
proporcionou um desempenho superior em termos de conforto térmico em comparação a 
outros arquivos climáticos analisados, com a zona permanecendo dentro dos parâmetros 
estabelecidos pela norma ISO 9241 por 2.960 horas ao longo do período anual, o que 
representa 33,8% do tempo. 

É fundamental destacar que a norma ISO 9241 estabelece parâmetros para o 
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conforto térmico em ambientes, considerando fatores como temperatura, umidade e 
velocidade do ar. Assim, o fato de a zona ter permanecido dentro desses parâmetros 
por um tempo significativo indica que o ambiente oferece condições adequadas para a 
realização de atividades. 
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RESUMO: Este trabalho apresenta as 
considerações sobre uma prática de 
campo realizada com estudantes do curso 
de Ciências Biológicas na Ilha Anchieta, 
região de Ubatuba/SP. O trabalho de 
campo consistiu em observação, descrição 
e interpretação da sucessão ecológica e 
das feições geomorfológicas e geológicas 
presentes ao longo das trilhas oficiais 
existentes na Ilha, além de sua evolução 
histórica. Dada a interdisciplinaridade da 
atividade desenvolvida, os alunos tiveram 
acesso a uma grande quantidade de 
informações, em situação e condições 
que não seriam possíveis de serem 
atingidas em sala de aula, significando 
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substancialmente o ensino-aprendizagem 
dos temas discutidos. 

PALAVRAS-CHAVE: Trabalho de campo. 
Ensino de Biologia. Ensino de Geociências. 
Ilha Anchieta. 


ABSTRACT: This paper presents the 
considerations about a field practice 
carried out with students of the Biological 
Sciences course in Anchieta Island, region 
of Ubatuba/SP. The fieldwork consisted of 
observation, description and interpretation 
of the ecological succession and the 
geomorphological and geological features 
present along the official trails existing on the 
Island, in addition to its historical evolution. 
Given the interdisciplinarity of the activity 
developed, the students had access to a 
large amount of information, in a situation 
and conditions that would not be possible to 
be reached in the classroom, substantially 
meaning the teaching-learning of the topics 
discussed. 

KEYWORDS: Fieldwork. Biology teaching. 
Teaching of Geosciences. Anchieta Island. 


INTRODUÇÃO 


O ensino atualmente praticado 
já não é mais o mesmo de outrora. Nas 
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últimas décadas, muitas discussões sobre como organizar o fazer pedagógico, de forma 
a torná-lo mais significativo aos alunos, tem norteado as pesquisas na área educacional e 
provocado reflexões sobre a prática docente (BELIZARIO, 2020). Nesse sentido, cresce a 
busca por metodologias ativas e práticas de ensino que promovam a interdisciplinaridade 
e o desenvolvimento do pensamento crítico ao estudante, aumentando seu protagonismo, 
seu interesse e seu envolvimento com a disciplina estudada. 

Uma prática comumente utilizada nos cursos de Ciências Biológicas, e nas Ciências 
Naturais em geral, é o trabalho de campo. Os trabalhos de campo são aquelas atividades 
que envolvem o deslocamento dos alunos para um ambiente não formal e alheio aos 
espaços de estudos contidos na escola (FERNANDES, 2007; SALVETTI e FRENEDOZO, 
2022). Esse tipo de atividade é adotado tanto no ensino superior quanto na educação 
básica, sendo utilizado como prática na produção ou reprodução de conhecimentos por 
professores, pois possibilitam a observação e contato direto com elementos da paisagem 
e processos que ocorrem no espaço geográfico (SALVETTI e FRENEDOZO, 2022), 
favorecendo a compreensão de fenômenos naturais (SANTOS & JACOBI, 2011) pela 
interação do estudante com situações reais que estimulam sua curiosidade e sentidos 
(VIVEIRO e DINIZ, 2009, SCHRADER e FRENEDOZO, 2015) e desenvolvendo um 
processo criativo relevante para a formação do estudante (CARNEIRO et al., 1993). 

Os trabalhos de campo como atividades de ensino-aprendizagem podem 
proporcionar aos alunos a motivação de descobrir e aprender os conteúdos curriculares 
a partir da avaliação de um espaço natural existente ou de sua própria realidade, externo 
ao ambiente escolar formal (GRANDI e MOTOKANE, 2012; SCHRADER e FRENEDOZO, 
2015). A aplicação de trabalhos de campo possibilita o ensino das temáticas físico-naturais 
e ambientais, por muitas vezes muito difíceis de serem exploradas em sala de aula, tanto 
pela limitação do livro didático quanto pela acomodação e a fragilidade teórica docente 
(BELIZARIO, 2020). 

Entendendo o trabalho de campo como uma metodologia ativa aplicada em um 
espaço não-formal de ensino, a atuação do professor torna-se essencial nesse processo, 
pois com sua abordagem possibilita aos alunos exercitarem sua autonomia, dando 
importância a atuação deles na construção dos seus próprios conhecimentos, valorizando 
suas expectativas e seus conhecimento prévios (SANTOS e JACOBI, 2011; SCHRADER e 
FRENEDOZO, 2015; BELIZARIO, 2020). Além disso, contribui para o desenvolvimento do 
pensamento crítico, pois faz com que os alunos exerçam a sua capacidade de posicionar- 
se, de elaborar e aplicar projetos, de trabalhar em grupo de maneira interativa e produtiva, 
além de relacionar o que se está aprendendo, o que se sabe e a prática (BARBOZA e 
RODRIGUES, 2016; BELIZARIO, 2020). 

Com isso em vista, o presente trabalho relata uma prática de campo realizada 
com estudantes do curso de Ciências Biológicas para estudo de conceitos relacionados 
a ecologia, geologia, dinâmica costeira e impactos ambientais de uma área litorânea 
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localizada no município de Ubatuba/SP. 


TRABALHOS DE CAMPO 


Trabalho de campo é toda atividade realizada externamente ao ambiente formal 
escolar, e que proporciona ao estudante o contato com múltiplas realidades (GRANDI 
e MOTOKANE, 2012). Tradicionalmente, os trabalhos de campo são aplicados nas 
disciplinas que desenvolvem estudos do meio (natural ou artificial), e se concretiza pela 
imersão orientada pelo professor do estudante na complexidade de um determinado espaço 
geográfico (LOPES e PONTUSCHKA, 2009), permitindo o estabelecimento de um diálogo 
inteligente com o mundo e a produção de conhecimento científico (GRANDI e MOTOKANE, 
2012). 

Durante uma atividade de campo, os alunos são convidados a explorar, por si 
próprios, a realidade do meio em que estão inseridos, analisando os fenômenos (naturais, 
ecológicos, sociais, urbanos etc.) que ocorrem no mundo real, e a comparar essas 
informações com os conceitos obtidos em sala de aula ou na literatura, alicerçando a 
consolidação do conhecimento e a aculturação científica. 

Nesse sentido, os trabalhos de campo se tornam um exemplo interessante 
de metodologia ativa de prática de ensino desenvolvida em um espaço não-formal 
de aprendizagem, que prima pela atividade tanto do professor quanto do aluno. As 
metodologias ativas diferem da “Pedagogia Tradicional” ao dar mais ênfase ao aluno em 
detrimento do professor, fazendo com que o professor deixe simplesmente de transmitir 
os conteúdos (prontos e acabados) e passe a “dirigir, estimular e orientar as condições 
externas e internas do ensino, de modo que, pela atividade dos alunos os conhecimentos 
e habilidades façam progredir seu desenvolvimento mental” (LIBÂNEO, 2013, p. 113). 
Fundamentam as metodologias ativas: a aprendizagem significativa, a interação do aluno 
com o outro (professor e demais alunos) e com o objeto de conhecimento, a valorização da 
experiência do indivíduo e a autonomia do sujeito (DIESEL, BALDEZ e MARTINS, 2017). 

Ao entrarmos em contato com o mundo real, as possibilidades de apreensão dos 
conhecimentos científicos se ampliam. COMPIANI (2007, p. 35) deixa claro que “em relação 
ao ensino [das Ciências Naturais], o campo — a natureza — significa o contato com objetos, 
fenômenos concretos, ambiente, ou seja, o contexto a partir do qual se criam situações e 
estratégias de aprendizagem”. Assim, a partir do contato direto com os fenômenos reais 
que compõem o meio natural, podemos programar atividades que explicitem o processo de 
formação desses fenômenos, permitindo a melhor compreensão da relação do ser humano 
com a natureza (COMPIANI, 1990) e das ações do Homem na construção das paisagens 
(MUNHOZ, 2013), contribuindo para que os alunos percebam sua real participação no 
mundo. 

Desta forma, o conhecimento científico adquirido pelas práticas de campo 
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passa a se relacionar diretamente à atuação da sociedade na natureza, uma vez que a 
compreensão da ciência por todos permite que o sujeito perceba sua participação real no 
mundo (COMPIANI; 1990; BONITO, 1999; GRANDI e MOTOKANE, 2012; MUNHOZ, 2013; 
LEMOS, 2021). 

Além disso, é inegável que as práticas de campo despertam nos estudantes a 
curiosidade, pois o simples fato de estarem fora do ambiente formal escolar já é considerado 
um incentivo para que eles participem ativamente da atividade proposta. 

Exatamente por isso, as atividades de estudos do meio possuem, pois, papéis 
didáticos importantíssimos (LOPES e PONTUSCHKA, 2009; LEMOS, 2021), e devem ser 
pautadas nas “funções que determinada atividade assume dentro do processo de ensino- 
aprendizagem, decididas de maneira deliberada ou não, que exercem algum significado 
para o alcance de objetivos didáticos” (COMPIANI e CARNEIRO, 1993, p. 90). Desse 
modo, as atividades não devem se resumir a deixar, momentaneamente, o espaço formal 
de aprendizagem, a sala de aula, e se transformar num mero passeio pelo espaço informal, 
mas precisam ter um significado para o aprendizado do aluno de modo a contribuir com 
suas leituras de mundo e de sua formação como cidadão crítico (MUNHOZ, 2013). 

Torna-se claro, portanto, que as práticas de campo devem ser bem planejadas e 
direcionadas, de modo a desenvolver nos alunos as habilidades e competências esperadas 
(CRUZ, 1997; LOPES e PONTUSCHKA, 2009; MUNHOZ, 2013; SILVA, FARIAS e LEITE, 
2019). A seleção do lugar a ser visitado, a determinação dos objetivos a serem alcançados 
e questões a serem respondidas na pesquisa de campo, o planejamento, a execução e 
a avaliação da atividade, devem ser orientadas tanto pela “dialogicidade” quanto pelo 
despertar da “curiosidade epistemológica” de todos os membros da comunidade escolar 
envolvida (FREIRE, 2000). Ou seja, todas as etapas e ações que estruturam uma prática 
de campo devem ser mediadas pelos professores e realizadas na busca de acordos e 
contratos pedagógicos que têm, como ponto de partida e chegada, a realidade do espaço 
geográfico estudado (LOPES e PONTUSCHKA, 2009). 

Posteriormente, é importante ainda avaliar as práticas aplicadas de modo a se 
determinar se os objetivos propostos com aquela atividade foram atingidos. 


TRABALHO DE CAMPO NO PARQUE ESTADUAL DA ILHA ANCHIETA, 
UBATUBA/SP 

Com o interesse de apresentar a alunos do curso superior de Ciências Biológicas, 
foi organizado um trabalho de campo no Parque Estadual da Ilha Anchieta, constituído por 
uma área insular de proteção ambiental administrada pelo Instituto Florestal do Governo do 
Estado de São Paulo, localizada em Ubatuba, litoral norte do Estado de São Paulo (figura 


1). 
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Google Earth 


Figura 1. Localização do Parque Estadual da Ilha Anchieta. 


Fonte: do autor. 


A Ilha Anchieta localiza-se na parte interna da plataforma continental e é separada 
do continente por um estreito canal chamado “Boqueirão”, de 600 m de largura e 35 m 
de profundidade. Todo o Parque Estadual está inserido no Domínio da Floresta Atlântica 
e abriga remanescentes de Floresta Ombrófila Densa e de Restinga, de clima tropical 
úmido (CICCHI et al., 2009). Como resultado de fortes pressões antrópicas, a Ilha Anchieta 
compreende hoje um mosaico vegetacional composto por florestas em diferentes estágios 
de regeneração e formações campestres antrópicas (CICCHI et al., 2009), o que resultou 
em perda expressiva de sua biodiversidade original. Geologicamente, a ilha compõe 
os Domínios Costeiros (no conceito de HASUI et al., 1981), representado por terrenos 
metamórficos migmatíticos com granitização variada semelhante aos observados na área 


continental. 
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O Parque possui algumas trilhas pré-definidas, cujo percurso é realizado sob a 
orientação de guias do local. Para esta prática foram percorridas duas trilhas principais. A 
primeira parte da Praia do Presídio, localizada na face norte da Ilha (figura 2) e termina na 
Praia do Sul, na porção ocidental da Ilha. A segunda trilha parte da Praia do Sul e vai até a 
Praia das Palmas, passando pelo Cemitério dos Búlgaros (figura 2). 


a: 
'aia das is 


Praia do Presídio 


Presídio da Ilha Anchieta 


Figura 2. Trilhas realizadas durante os trabalhos de campo na Ilha Anchieta. 


Fonte: do autor. 


Como a Ilha foi fortemente afetada pela ação antrópica ao longo do século XX, 
sobretudo devido ao funcionamento do presídio, a vegetação da ilha apresenta reduzida 
diversidade florística (GUILLAUMON e FONTES,1992) e a fauna nativa foi praticamente 
eliminada, sendo boa parte das espécies ali presentes introduzidas em 1983 pelo Zoológico 
de São Paulo (GUILLAUMON et al. 1989), resultando em grave desequilíbrio ecológico. 

Nesse contexto ecológico, os alunos puderam observar, ao longo das trilhas 
percorridas, a sucessão vegetal e identificar algumas espécies da flora e fauna existentes. 
Muitas das espécies vegetais identificadas foram consideradas exóticas para o local, pois 
não pertencem a mata nativa da Ilha ou mesmo à Mata Atlântica. 

Na Praia do Sul, os alunos participaram de mergulhos guiados para observação da 
fauna presente no costão rochoso submerso (figura 3-1) e aéreo, além da faixa de areia da 
praia (figura 3-3). Nesse local foi possível observar, ainda, as características da dinâmica 
costeira e fluvial da região (figura 3-4). 


Estudos multidisciplinares: Ciências exatas e da terra e engenharias Capítulo 5 


Figura 3. Atividades realizadas na Praia do Sul, Ilha Anchieta. 1. Alunos participando de mergulho 
guiado para observação do costão rochoso submerso. 2. Crustáceo filtrador (Chthamalus) fixos nas 
rochas do costão rochoso. 3. Caranguejo (Ocypode quadrata - Maria-farinha ou Guaruça) identificado 
na faixa de areia da Praia do Sul. 4. Canal fluvial desaguando na Praia do Sul. 


Fonte: do autor. 


Atrilha 2, que corta a área de restinga próxima a Praia das Palmas, é toda construída 
sobre plataformas suspensas (figura 4-1), o que permite a visão detalhada da vegetação 
e dos animais presentes nesse trecho. Essa trilha (vide figura 2) passa pelo Cemitério dos 
Búlgaros (figura 4-2), importante registro arqueológico existente na Ilha. 
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Figura 4. Trecho da Trilha 2, sobre plataformas suspensa de madeira. 1. Alunos observam explicação 
do guia sobre a vegetação presente no local. 2. Vista geral do Cemitério dos Búlgaros. 


Fonte: do autor. 


O término da visita guiada se deu na área do Presídio da Ilha Anchieta (figura 5), 
local histórico que data do início do século XX e que foi utilizado como presídio político e de 
segurança máxima até seu fechamento, em 1955 (MEMORIAL DA RESISTÊNCIA, 2022). 


Figura 5. Vista geral da área do Presídio da Ilha Anchieta. 


Fonte: do autor. 


AVALIAÇÃO DA ATIVIDADE 


Ao término das atividades de campo os alunos foram convidados a dar sua opinião 
sobre a prática realizada. Através de um questionário tipo likert, onde os alunos indicavam 
se concordavam ou discordavam das afirmações propostas, foram apresentadas as 
seguintes afirmativas: 

1. Você ficou satisfeito com o trabalho de campo realizado na disciplina. 


2. A atividade de campo executada contribuiu para consolidar os conceitos teóricos 
trabalhados em sala de aula. 


3. A aula de campo foi adequada e seguiu os roteiros propostos. 
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4. Aprendi novas técnicas de levantamento de dados durante a aula de campo. 


5. Antes da aula de campo, não me considerava apto(a) para executar as tarefas 
solicitadas. 


6. Após a aula de campo, me considero mais preparado para executar levantamentos 
de dados na prática. 


7. Os professores e monitores que conduziram a aula de campo possuem 
conhecimento sobre os tópicos discutidos. 


8. Atividades de campo são práticas importantes para a minha formação profissional? 


As respostas obtidas estão sintetizadas na figura 6 a seguir. 

Em termos gerais, é possível notar que a maioria dos alunos que responderam ao 
formulário concordam que a atividade de campo realizada foi importante para consolidar os 
conceitos e informações discutidos em sala de aula, e contribuir para a formação profissional 
desses estudantes. Adicionalmente, os alunos afirmaram que o trabalho de campo permitiu 
maior contato com o mundo real, possibilitando visualizar os conceitos trabalhados em sala 
de aula fora do ambiente escolar, considerado às vezes deveras controlado e maçante. 


5 - concordo totalmente 


iii mm 
FTTT—=" ———— 


3 - não estou decidido 
2 - discordo 


1 - discordo totalmente 


(0) 1 2 3 4 5 6 A 8 9 
mAtividades de campo são práticas importantes para a minha formação profissional 
EH Os professores e monitores que conduziram a aula de campo possuem conhecimento sobre os tópicos discutidos 
E Após a aula de campo, me considero mais preparado para executar levantamentos de dados na prática 
E Antes da aula de campo, não me considerava apto(a) para executar as tarefas solicitadas 
E Aprendi novas técnicas de levantamento de dados durante a aula de campo 
EA aula de campo foi adequada e seguiu os roteiros propostos 
EA atividade de campo executada contribuiu para consolidar os conceitos teóricos trabalhados em sala de aula 


E Fiquei satisfeito com o trabalho de campo realizado na disciplina 


Figura 6. Respostas obtidas dos alunos sobre a prática efetuada. 


Fonte: do autor. 
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Dada a interdisciplinaridade da prática efetuada, considera-se que a atividade 
desenvolvida permitiu aos alunos o acesso a uma grande quantidade de informações, 
em situação e condições que não seriam possíveis de serem atingidas em sala de aula, 
significando substancialmente o ensino-aprendizagem dos temas discutidos. 


CONCLUSÃO 


O emprego dos trabalhos de campo como metodologia ativa de ensino se mostra 
uma ferramenta muito eficaz para a consolidação do conhecimento dentre os estudantes, 
por permitir que eles, investigando e avaliando o meio natural, avaliem suas feições e 
características e cheguem a conclusões sobre o ambiente estudado. 

As atividades de campo permitem que os estudantes tenham total acesso às feições 
do meio natural, situação impossível dentro do ambiente escolar. Assim, entende-se que 
os trabalhos de campo são a expressão máxima do emprego dos ambientes não-formais 
na aprendizagem e possuem enorme potencial que pode ser explorado para consolidar os 
conceitos trabalhados pelos professores em sala de aula. 

Complementarmente, esse tipo de prática desperta no aluno o pensamento e o senso 
crítico, despertando o interesse pelos estudos ambientais e estimulando-os a desenvolver 
técnicas de levantamento de campo mais aprimoradas. 

Por fim, os próprios estudantes sentem-se interessados em participar das atividades 
de campo, e julgam esse tipo de prática importante para sua vivência acadêmica e 
consolidação do conhecimento científico. 
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